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Introduction 


Amendment 22 to the Gulf of Mex- 
ico Fishery Management Council’s 
(GMFMC) Reef Fish Fishery Manage- 
ment Plan (GMFMC') dictates manda- 
tory observer coverage. In July 2006, 
in collaboration with the commercial 
fishing industry and the GMFMC, 
the National Marine Fisheries Ser- 
vice’s (NMFS) Southeast Fisheries 
Science Center (SEFSC) implemented 
a mandatory observer program to char- 
acterize the commercial reef fishery 
operating in the U.S. Gulf of Mexico 
(Gulf). 


'GMFMC. 2005. Amendment 22 to the Reef Fish 
Management Plan. Gulf Mex. Fish. Manage. 
Counc., Tampa, Fla. (available at http://www. 
gulfcouncil.org). 


The authors are with the Southeast Fisheries Sci- 
ence Center, National Marine Fisheries Service, 
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Descriptions of the U.S. Gulf of Mexico Reef Fish 
Bottom Longline and Vertical Line Fisheries Based on Observer Data 
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JEFF R. PULVER, REBECCA C. SMITH, and JO ANNE WILLIAMS 


This fishery consists of approxi- 
mately 890 Federally permitted vessels 
(SERO?). Primary gears used include 
bottom longline, vertical line (bandit or 
handline), and more recently, modified 
buoy gear. Although many reef fish spe- 
cies are retained, the predominant target 
species are groupers, Epinephelus spp.. 
and snappers, Lutjanus spp. 
Longliners off the coast of Florida 
generally target red grouper, Epineph- 
elus morio, in shallow waters, and in 
deeper waters yellowedge grouper, E. 
flavolimbatus; tilefish (Malacanthidae), 
and sharks (Carcharhinidae). Vertical 
line vessel operators target shallow- 
water grouper (e.g. red grouper), red 
snapper, Lutjanus campechanus, and 
may also seek yellowedge grouper and 
vermilion snapper, Rhomboplites au- 
rorubens. From historical effort data, 


*SERO. 2010. Fishery permits and fishery 
quotas. Southeast Reg. Off., Natl. Mar. Fish. 
Serv., NOAA, St. Petersburg, Fla. (available at 
http://sero.nmfs.noaa.gov). 


ABSTRACT—In July 2006, a manda- 
tory observer program was implemented to 
characterize the commercial reef fish fish- 
ery operating in the U.S. Gulf of Mexico. 
The primary gear types assessed included 
bottom longline and vertical line (bandit 
and handline). A total of 73,205 fish (183 
taxa) were observed in the longline fishery. 
Most (66%) were red grouper, Epineph- 
elus morio, and yellowedge grouper, E. 
flavolimbatus. Ji the vertical line fishery, 
89,015 fish (178 taxa) were observed of 
which most (60%) were red snapper, Lut- 
janus campechanus, and vermilion snapper, 
Rhomboplites aurorubens. Based on sur- 
face observations of discarded under-sized 
target and unwanted species, the major- 
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ity of fish were released alive; minimum 
assumed mortality was 23% for the vertical 
line and 24% for the bottom longline fish- 
ery. Of the individuals released alive in the 
longline fishery, 42% had visual signs of 
barotrauma stress (air bladder expansion/ 
and or eyes protruding). In the vertical line 


fishery, 35% of the fish were released in a 


stressed state. Red grouper and red snap- 
per size composition by depth and gear type 
were determined. Catch-per-unit-effort for 
dominant species in both fisheries, illus- 
trated spatial differences in distribution 
between the eastern and western Gulf. 
Hot Spot Analyses for red grouper and red 
snapper identified areas with significant 
clustering of high or low CPUE values. 


most commercial fishing effort for red 
snapper occurs in the western Gulf of 
Mexico (SEDAR?). 

In November 1984, the Reef Fish 
Fishery Management Plan (GMFMC‘*) 
was implemented to rebuild declining 
reef fish stocks. Since that time, Federal 
regulations have restricted size and 
landings of several reef fish species. 
Weight quotas regulate commercial 
landings for grouper, with 7.57 mil- 
lion Ibs for shallow-water grouper and 
1.02 million lbs for deepwater grouper 
(SERO?). The current total allowable 
catch (TAC) for red snapper is 6.3 mil- 
lion Ibs, divided between the commer- 
cial (51%) and recreational (49%) fish- 
ing sectors. An individual fishing quota 
(IFQ) program for the commercial red 
snapper fishery was implemented in 
2007 and for the grouper and tilefish 
fisheries in 2010. 

Certain areas for reef fish are 
closed or restricted based on gear type 
(GMFMC>). Federal waters are closed 
in the Tortugas North and Tortugas 
South Ecological Reserves in the Florida 
Keys National Marine Sanctuary and the 
Madison and Swanson and Steamboat 
Lumps Marine Reserves off the west 
central Florida coast. Longline and 
other buoy gear are prohibited inside 


SSEDAR. 2005. Stock assessment report of 
SEDAR 7 Gulf of Mexico red snapper. Southeast 
Data Assessment and Review, South Atl. Fish. 
Manage. Counc., Charleston, SC (available at 
www.sefsc.noaa.gov/sedar/). 

4GMFMC. 1984. Reef Fish Management Plan. 
Gulf Mex. Fish. Manage. Counc., Tampa, Fla. 
(available at http://www.gulfcouncil.org). 
>GMFMC. 2010. Commercial fishing regulations 
for Gulf of Mexico Federal waters. Gulf Mex. 
Fish. Manage. Counc., Tampa, Fla. (available at 
http://www.gulfcouncil.org). 
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the 50-fm contour west and the 20-fm 
contour east of Cape San Blas, Fla. 

In May 2009, an emergency rule to 
protect sea turtles (Cheloniidae and 
Dermochelyidae) went into effect pro- 
hibiting the use of bottom longline gear 
east of Cape San Blas, Fla., shoreward of 
the 50-fm contour. Modification through 
subsequent regulations (GMFMC>) 
prohibited bottom longline gear east 
of Cape San Blas, Fla., shoreward of 
the 35-fm contour from June through 
August, restricted the number of hooks 
onboard to 1,000, of which only 750 
could be rigged for fishing, and reduced 
the number of vessels through an en- 
dorsement system based on documen- 
tation of an average annual landing of 
at least 40,000 Ibs during 1999 through 
2007. 

The effectiveness of quota systems, 
size limits, and area closures as manage- 
ment tools has been debated (Coleman et 
al., 2000; Nieland et al., 2007; Stephen 
and Harris, 2010). Once a vessel’s red 
snapper quota is reached, for example, 
the vessel simply targets other reef fish, 
making red snapper a bycatch species. 
Currently, the minimum legal size for 
red snapper is 13 in total length (TL). 
The minimum size limit for red grouper 
was reduced from 20 in TL to 18 in TL, 
effective 18 May 2009 (GMFMC°). 

The mortality rates of both undersized 
target species and nontargeted species 
caught on the various gear types remains 
a pressing concern. Findings from mark- 
release mortality studies (Gitschlag and 
Renaud, 1994; Schirripa and Legault®; 
Burns et al.’) indicate variable rates of 
mortality based on depth and method 
of capture. 

In December 1993, SEFSC’s Galves- 
ton Laboratory implemented a volun- 
tary observer program to characterize 


*Schirripa, M. J., and C. M. Legault. 1999. 
Status of red snapper in U.S. waters of the 
Gulf of Mexico: updated through 1998. U.S. 
Dep. Commer., NOAA, Natl. Mar. Fish. Serv., 
Panama City Lab. Sustainable Fish. Div. Contrib. 
SFD-99/00-75. 

7Burns, K. M., N. F. Parnell, and R. R. Wilson, Jr. 
2004. Partitioning release mortality in the under- 
sized bycatch: Comparison of depth vs. hooking 
effects. MARFIN Grant No. NA97FF0349, 36 
p., on file at Southeast Reg. Off., Natl. Mar. Fish. 
Serv., NOAA, St. Petersburg, Fla. 


the fish trap, bottom longline, and 
bandit reel fisheries in the U.S. Gulf 
of Mexico (Scott-Denton and Harper’®; 
Scott-Denton’). Observer coverage 
of the commercial reef fish fishery 
operating primarily off the west coast 
of Florida and, to a lesser extent, off 
Louisiana, was conducted from 1993 
through 1995. Data from 576 sets 
aboard fish trap vessels, 317 sets from 
bottom longline, and 580 sets from 
bandit reel vessels were analyzed. 
Findings from this study revealed a 
low proportion (<5% of total number 
caught) of fish discarded dead (im- 
mediate mortality) based on surface 
observations. However, due to the 
number of fish released in stressed 
state (air bladder expansion and/or eyes 
protruding), total predicted red snapper 
discards of 25% to 30% were used to 
estimate the number of discarded fish 
at age that died and thus contributed to 
fishing mortality (Goodyear'®). 

The continuing goal of the current 
observer program is to provide quan- 
titative biological, vessel, and gear- 
selectivity information relative to the 
directed reef fish fishery. The specific 
objectives are to: 1) provide general 
fishery bycatch characterization for 
finfish species taken by this fishery, 2) 
estimate managed finfish discard and 
release mortality levels, and 3) esti- 
mate protected species bycatch levels. 
The specific objectives of this report 
are to: 1) summarize trip, vessel, en- 
vironmental, and gear characteristics, 
2) quantify fish and protected species 
composition and disposition based on 
surface observations, 3) examine size 
composition of target species, and 4) 
estimate catch-per-unit-effort (CPUE) 


‘Scott-Denton, E., and D. Harper. 1995. Char- 
acterization of the reef fish fishery of the eastern 
Gulf of Mexico. SEFSC Rep. to Gulf Fish. Man- 
age. Counc. July 17, 1995, Key West, Fla., 45 p. 
°*Scott-Denton, E. 1996. Characterization of 
the reef fish fishery of the eastern U.S. Gulf of 
Mexico. MARFIN Grant No. 95MFIHO7. Suppl. 
Rep. to MARFIN Grant No. 94MARFINI7, on 
file at Southeast Reg. Off., Natl. Mar. Fish. Serv., 
NOAA, St. Petersburg, Fla. 

Goodyear, C. P. 1995. Red snapper in U.S. 
waters of the Gulf of Mexico. U.S. Dep. 
Commer., NOAA. Natl. Mar. Fish. Serv., South- 
east Fish. Sci. Cent., Miami Lab. Rep. Contrib. 
MIA 95/96-05, 171 p. 


trends and spatial distribution for domi- 
nant species. 


Methods 


Protocol sampling modification, ran- 
domized vessel selection, and observer 
deployment through mandatory efforts 
began in 2006 for the commercial reef 
fish fishery. NMFS observers were 
placed on reef fish vessels operating 
throughout the Gulf of Mexico based 
on randomized selection stratified by 
season, gear, and region. Proportional 
sampling effort, based on coastal log- 
book data, among seasons and gears 
in the eastern and western Gulf of 
Mexico was recommended by SEFSC 
stock assessment scientists in 2006 
and used thereafter for vessel selection 
stratification purposes using annual 
updated effort data. Thus, proportional 
sampling was used to direct coverage 
levels (based on sea days, the National 
metric for percent observer coverage 
levels) toward region and gear strata 
with higher levels of fishing effort, while 
continuing to sample strata with lower 
fishing effort. 

In 2008, for the longline fishery, 
seven trips were not selected through 
the mandatory process. Instead the trips 
were based on voluntary cooperation 
as part of a pilot project to assess the 
effectiveness of electronic monitoring 
equipment. Observers placed on these 
vessels were equipped with closed- 
circuit video cameras and associated 
electronics. Results of this study are 
reported by Pria et al. (2008). 

In February 2009, increased coverage 
was directed toward the bottom longline 
fishery in the eastern Gulf to monitor 
for sea turtle interactions. In response 
to the bottom longline closure inside 
the 50-fm contour in the eastern Gulf in 
2009, some traditional longline vessels 
used modified buoy gear. This gear type 
was deployed during three trips inside 
50 fm in December 2009 with observ- 
ers onboard. 

Shrimp statistical zones (Patella, 
1975) were used to delineate area 
designations (Fig. 1). Conventionally, 
statistical areas 1-9 represent areas off 
the west coast of Florida, 10-12 delin- 
eate Alabama/Mississippi, 13—17 depict 
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Louisiana, and 18-21 denote Texas. For 
the reef fish fishery, statistical areas 1-8 
represent the eastern Gulf and areas 
9-21 the western Gulf. Seasonal catego- 
ries were: January through March, April 
through June, July through September, 
and October through December. The 
three primary gear types assessed in- 
cluded bottom longline, bandit reel, and 
handline. The latter two were combined 
to represent the vertical line fishery. 
Among the several provisions pro- 
mulgated under Magnuson-Stevens 
Conservation and Management Act 
(MSFCMA) § 303(b)(8) is the man- 
date for Federal permit holders to have 
a current Commercial Fishing Vessel 
Safety Examination decal prior to the 
selection period for mandatory observer 
coverage. The safety decal requirement, 
in combination with other factors, led 
to low vessel compliance, especially in 
the first 2 years of the study. A dedicated 
effort by NOAA Office of Law Enforce- 
ment (OLE) has substantially increased 
compliance (>95%). Additionally, a 
minimum sea day requirement by gear 
type was established to prevent early 
trip termination due to observer effect. 
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Figure |.—Distribution of sampling effort (sets) based on observer coverage of the U.S. Gulf of Mexico bottom longline reef 
fish fishery from August 2006 through November 2009. 


Reef fish permit holders are required 
to carry an observer for a minimum of 
7 days during a selection period when 
using longline gear, 3 days for bandit 
gear, and 2 days for handline. 

Once deployed, vessel length, hull 
construction material, gross tonnage, 
engine horsepower, and crew size were 
obtained for each vessel. For each set 
(the location of gear placement at a 
defined time), the type, number, and 
construction material of the fishing 
gear were recorded. Latitude, longitude, 
depth, and environmental parameters 
including sea state and bottom type 
were recorded at the start of each set. 
The total time the gear remained in 
the water (soak or fishing time) was 
calculated. 

Fishery data were obtained from each 
set. If a set could not be sampled due 
to time constraints or weather condi- 
tions, a minimum of location, depth, 
and fishing time were recorded. The 
condition of fish when brought onboard 
was categorized into one of the follow- 
ing: 1) live—normal appearance, 2) 
live—stomach/air bladder protruding, 
3) live—eyes protruding, 4) live—com- 
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bination of 2 and 3, 5) dead on arrival, 
or 9) not determined.'! Categories 2 
through 4 were combined to represent 
a stressed condition. 

Fate of fish after release was recorded 
as alive if it swam down or as discarded 
dead if it swam erratically, floated, or 
sank, or if undetermined. Nontarget 
and undersized target species were pro- 
cessed first by recording length, weight, 
condition when brought onboard, and 
fate after release to provide an estimate 
of immediate mortality (number dis- 
carded dead divided by the number of 
total discards). 

If venting occurred, air bladders of 
live discarded fish were punctured in 
the same manner as demonstrated by the 
captain and crew if requested. Retained 
species were processed by recording 
length, weight, condition when brought 
onboard, and if kept or retained for bait. 
Sightings or captures of sea turtles were 
recorded in accordance with SEFSC 
protocol (NMFS, 2008). Data pertaining 
to sea turtle interactions were reported 


"Category 9 is the default for a condition that is 
unknown or not recorded. 
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to SEFSC for annual sea turtle mortality 
estimates. 

On some (19%) vertical line sets, due 
primarily to time constraints and the 
magnitude of the catch, not all reels were 
sampled for the set. The species total 
number was extrapolated proportionally 
based on subsampled reels for that set. 
Negative sets, or sets where no fish were 
caught, were included in CPUE calcu- 
lations. No extrapolation procedures 
were required for longline and modi- 
fied buoy sets (i.e. all hooks sampled). 

Overall catch rates are presented 
collectively for all years, areas, sea- 
sons, and depths. Due to data confi- 
dentiality rules, a minimum of three 
vessels were required for spatial and 
temporal stratification purposes, and 
analysis of modified buoy gear data 
was restricted. 

Effort was calculated using methods 
described by McCarthy and Cass- 
Calay.'* The number of hooks set at 
each location was multiplied by soak 
time to derive hook-hours. Catch rates 
were calculated in number of fish per 
hook-hour. For the vertical line fishery, 
total soak time was used for one set loca- 
tion using the sum of all hooks per reel. 
Therefore, effort may be overestimated 
due to the repeated deployment (e.g. 
drops) of multiple gear configurations 
(e.g. hooks) on the same reel at one set 
location. Moreover, average haul in time 
was not documented for all sets, there- 
fore not used in the effort calculation. 
For sets when the average haul in time 
was recorded, the average value was less 
than one minute. 

Ratio estimation was used for analy- 
ses of species-specific catch rates. As 
described by Snedecor and Cochran 
(1967) and Watson et al. (1999), the 
ratio estimation (1) below was used as 
the sample estimate of the mean. 


McCarthy, K. J., and S. Cass-Calay. 2006. 
Standardized catch rates for red grouper from the 
United States Gulf of Mexico handline, longline, 
and trap fisheries, 1990-2005. SEDAR 12-DW- 
16. Southeast Data Assessment and Review, 
South Atl. Fish. Manage. Counc., Charleston, SC 
(available at www.sefsc.noaa.gov/sedar/). 


where: R = ratio estimate, 

Y= extrapolated number for 
species of a particular dis- 
position code for selected 
strata, and 

X = hook-hours for selected 
strata. 


The estimated standard error of the 
estimate is given in equation 2: 


Ry = (2) 
x\ 


n(n-1) 
where: x= mean of hook-hours for 
selected strata, and 
number of sets occurring in 
selected strata. 


A density surface of CPUE, based 
on number of fish kept per 1,000 hook- 
hours for dominant species by fishery, 
was created using Fishery Analyst.!3-!4 
This is an ArcGIS extension developed 
to graphically present temporal and 
spatial trends in fishery statistics (Riolo, 
2006). A search radius of 25 km was 
used to ensure the search parameter 
encompassed the maximum length of a 
fishing set. A cell size of 5 km produced 
the desired resolution. 

Density of catch and effort values for 
each 5 km cell were calculated by sum- 
ming those values contained within the 
25 km search radius and dividing the 
value by the area of the circle as defined 
by the search radius. A summary CPUE 
value for all years combined was calcu- 
lated for each cell by calculating CPUE 
values for individual years and dividing 
by the number of years for which fishing 
activity occurred in that cell. 

To identify patterns in CPUE for the 
most frequently captured species in each 
fishery, a local spatial statistic, the Getis- 
Ord Gi* (Gi*), was calculated using the 
Hot Spot Analysis tool in ArcGIS!>, to 


'5Fishery Analyst, Mappamondo GIS, Via 
Rubens 3, 43100 Parma(PR)-Italy. 

'4Mention of trade names or commercial firms 
does not imply endorsement by the National 
Marine Fisheries Service, NOAA. 

'SArcGIS 9.3 Computer Software. 380 New 
York Street, Redlands, Calif. 92373. 


locate clusters of features with similarly 
high or low values. The Gi* statistic was 
also calculated for all discarded and kept 
species in order to assess if geographi- 
cal areas of particularly high levels of 
bycatch occurred. 

The Hot Spot Analysis tool evalu- 
ates each feature within the context of 
neighboring features. If the value of the 
feature is high, and the values for all of 
its neighboring features are also high, 
it is a part of a hot spot. Conversely, 
if a feature is surrounded by similarly 
low values, it is identified as a cold 
spot. The Gi* statistic is a Z-score test 
Statistic. For statistically significant 
positive Z-scores, the larger the Z-score 
is, the more intense the clustering of 
high values. The Z-score can produce 
misleading results when used with 
local statistics because the test assumes 
independence between features. Since 
the GIS runs the test to calculate a Z- 
score for each feature, the test will end 
up using many of the same neighbors 
for adjacent features (Mitchell, 2005). 
For this reason, the statistical tests as- 
sociated with local measures of spatial 
autocorrelation for data exploration 
were used, rather than as confirmatory 
statistical testing (Nelson and Boots, 
2008). 

To standardize bycatch (discard) 
estimates as prescribed in “Evaluating 
Bycatch” (NMFS, 2004), the coef- 
ficient of variation (CV) was used as 
a measure of precision for bycatch es- 
timates. CV estimates were calculated 
by dividing the estimated standard error 
by the estimate of the mean CPUE 
(number per hook-hour) for Federally 
managed discarded species. Less than 
0.3% of the total fish processed had an 
undetermined fate code and were as- 
sumed to be discarded in an unknown 
condition. 

Length data are given for the domi- 
nant target species. Fish measurements 
were recorded in metric units for age 
and growth assessment. To be consistent 
with the current regulatory mandates 
relative to size limits, metric measure- 
ments were converted to U.S. system 
equivalents. Fork to total length con- 
versions for red grouper were based on 
metric regression (Lombardi-Carlson 
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et al.'°). Red snapper total lengths were 
derived from fork length measurements 
using equation 3 (SEDAR, 2005): 


TL. 
0.1729 + FL (in) * 1.059. (3) 


After converting, length values were 
placed into | in intervals. Any lengths 
ranging from 19.000 to 19.999, for ex- 
ample, were categorized as 19 in. Hence, 
some degree of error is assumed. Only 
length measurements were considered. 
Weight data were not recorded for all 
specimens, therefore were not included 
in the analysis. 


Results 


Fishing Characteristics 


From July 2006 through December 
2009, data from 9,468 sets collected 
during 308 trips (1.919 sea days) aboard 
205 reef fish vessels were analyzed. 
Number of trips, sets, sea days, and 
percent coverage levels are given by 
year and project (Table 1). 

Trip, vessel, set, and gear character- 
istics varied by gear type (Tables 2, 3). 
Trip length averaged 11.7 days for long- 
line and 4.8 days for vertical line. Vessel 
length ranged from 23 to 70 ft, with 
longline vessels typically at the larger 
end of the range. The majority (=>85%) 
of vessels were fiberglass construction. 

For longline, the distance of mainline 
set at a location averaged 5.6 nmi. Mean 
gangion length was 5.8 ft. On average, 
991 circle hooks were set at a location. 
Most hooks (43%) were 13 aught in size 
and ranged from 12 to 15 aught. In the 
vertical line sector, the number of reels 
used at a set averaged 3.3. The majority 
(51%) of reels were electric. The number 
of hooks deployed during a set averaged 
26 hooks, with circle hooks deployed 
most often. The majority (43%) of 
hooks were smaller hooks (8 aught) as 
compared to longline. 


'6Lombardi-Carlson, L. A., G. R. Fitzhugh, and 
J. J. Mikulas. 2002. Red grouper (Epinephelus 
morio) age-length structure and description of 
growth from the eastern Gulf of Mexico: 1992- 
2001. U.S. Dep. Commer., NOAA. Natl. Mar. 
Fish. Serv., Southeast Fish. Sci. Cent., Contrib. 
Ser. 2002-06, 42 p. 


Table 1.—Reef fish trips, sets, and sea days by year and project from July 2006 to December 2009. 


Trips by Year and Project 


Electronic Buoy 

Year Bandit Handline Longline Monitoring Gear Total 

2006 30 8 12 50 

2007 72 25 11 108 

2008 34 19 5 65 

2009 28 21 33 3 85 

Total 164 73 61 3 308 

Sets by Year and Project 
Electronic Buoy 

Year Bandit Handline Longline Monitoring Gear Total 

2006 1,078 62 201 1,341 

2007 2,424 505 194 3,123 

2008 1,353 298 110 245 2,006 

2009 1,361 310 753 574 2,998 

Total 6,216 1,175 1,258 245 574 9,468 

Sea Days by Year and Project 
Electronic Buoy Industry Percent 

Year Bandit Handline Longline Monitoring Gear Total Sea Days Coverage 
2006 184 12 113 309 21,379 1.4 
2007 396 69 120 585 38,200 1.5 
2008 219 38 45 108 410 37,348 1.1 
2009 162 36 397 20 615 36,818 1.6 
Total 961 155 675 108 20 1,919 133,745 1.4 


Fishing and environmental condi- 


tions differed by gear type (Tables 2, 


3). Average fishing depth for longline 
sets was 51.5 fm. Fishing depths were 
shallower (27.3 fm) for vertical line. 
Average soak time was 5.1 h for long- 
line and 0.7 h for vertical line. Most 
sets (>47%) occurred over rock bottom 
in seas <2 ft during daylight hours for 
both gear types. 


Bottom Longline 
Allocation of Sampling Effort 


Data from 68 trips aboard 48 bottom 
longline vessels from August 2006 
through November 2009 were analyzed. 
The capture of 73,205 fish (Table 4) 
occurred during 1,503 sets deploying 
traditional longline gear (Fig. 1). For 
longline, | 431 sets had associated effort 
data (7,232 h; 1,395,320 hooks). Ap- 
proximately 90% of fishing effort, based 
on hook-hours, occurred in the eastern 
Gulf. The greatest concentration of effort 
(hook-hours) occurred in statistical areas 
3 through 5 (Fig. 2), with most (35%) in 
area 4. By season, 20% of the sets oc- 
curred from January through March; 52% 
April through June; 16% July through 
September; and 12% October through 
December for all years combined. 


Species Composition 


Of the 73,205 fish (183 taxa) caught 
on longline gear, 46% of the individu- 
als were kept, 35% were released alive, 
12% were discarded dead, 4% were 
discarded with an unknown condition, 
and 3% were retained for bait (Tables 5 
and 6). By number, red grouper domi- 
nated the catch composition at 56%. 
Yellowedge grouper comprised 10% of 
the catch, followed by blueline tilefish, 
Caulolatilus microps, at 5%; red snap- 
per, tilefish, Lopholatilus chamaeleon- 
ticeps, and Atlantic sharpnose shark, 
Rhizoprionodon terraenovae, each at 
3%. All other species combined consti- 
tuted 20% of the catch. 

By category, red grouper, yellow- 
edge grouper, tilefish, and blueline 
tilefish comprised the majority (82%) 
of the 33,335 individuals kept by 
longliners. Four species (red grouper, 
Atlantic sharpnose shark, smooth dog- 
fish, Mustelus canis; and red snapper) 
accounted for 83% of the released alive 
category. Of the 25.471 individuals 
released alive, 42% exhibited visual 
signs of stress, while 46% exhibited 
a normal appearance. Of the 2,414 
individuals used for bait, the species 
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caught and used most often for bait and red snapper comprised the major- Atlantic sharpnose shark (34%), and 
were king snake eel, Ophichthus rex _ ity (81%) of 9,037 individuals in the _ red snapper (27%). The fate of 2,948 
(29%), and palespotted eel, Ophichthus discarded dead category. Approximate individuals was undetermined. Of 
puncticeps (11%). Red grouper, blue- | minimum assumed mortality was: red these, approximately 77% were red 
line tilefish, Atlantic sharpnose shark, grouper (20%), blueline tilefish (76%), grouper. 


Table 2.—Trip, vessel, set, gear, and environmental characteristics observed in the longline fishery from August 2006 to November 2009. 


Longline 


Trip Vessel Set Gear Environmental 


783 Sea Days Length Soak time Mainline material: Water Depth: 
68 trips aboard 48 vessels Avg: 48.3 ft Avg: 5.1 h Cable (92%) Avg: 51.5 fathoms 
1,503 sets Range: 35 to 69 ft (+ 2.9 s.d.) Monofilament (7%) (+ 37.8 s.d.) 
(+8.4s.d.) Range: 0.9 to 32.2 h Test Eastern: 44.5 
Avg: 1,472.8 Ibs Western: 51.5 
(+ 784 s.d.) Range: 19.3 to 212.0 


Range: 310 to 4,000 Ibs 


Trip Length Hull Construction Mainline Gangion: Sea State: 
Avg: 11.7 days Fiberglass: 85% Avg length: 5.6 nmi Monofilament (99.9%) «Nylon (0.1%) 0 to 2 foot seas: 46% 
(+3.8s.d.) Steel: 10% (+ 2.0 s.d.) Avg length: 5.8 ft 3 to 5 foot seas: 35% 
Range: 4 to 20 days Fiberglass/wood: 4% Range: 0.9 to 12.0 nmi (+ 2.1 s.d) 6 to 8 foot seas: 17% 


Range: 2.5 to 11.0 ft 8+ foot seas: 2% 


Crew size Engine Horsepower Hooks: Bottom type: 
1 to 3 individuals (excluding captain) Avg: 277.1 hp Avg: 991.1 hooks Rock: 47% 
(+ 205.3 s.d.) (+ 426.4 s.d.) Unknown: 14% 
Range: 76 to 1400 hp Range: 150 to 2,500 hooks Shell: 16% 
Type: Circle hooks (100%), offset Coral: 10% 
(63.4%), straight (36.6%) Mud: 8% 
Shaft length avg 2.1 in Sand: 2% 
Distance between hooks: Boulder, clay, and grass: 
Avg: 22.5 ft (+ 13.0 s.d.) 1% each 


Range: 7.0 to 75.0 ft 
Size: 13 aught (43%) 
Range: 12 to 15 aught 
Brand: Mustad®: 82% 
Eagle Claw®: 18% 


Table 3.—Trip, vessel, set, gear, and environmental characteristics observed in the vertical line fishery from July 2006 to December 2009. 


Vertical Line 


Trip Vessel Set Gear Environmental! 


1,116 Sea Days Length Soak time Reel type: Water Depth: 
237 trips aboard 157 vessels Avg: 39.2 ft Avg: 0.7 hrs (+ 1.1 s.d.) Electric: 51.4% Avg: 27.3 fathoms (+ 15.8 s.d.) 
7,391 sets Range: 23 to 70 ft (+ 9.6 s.d.) Range: 0.02 to 15.3 h Hydraulic: 21.7% Range: 0.7 to 305.0 
Haul in time Hand: 27.0% 
Recorded: 68% 
Avg: 0.8 min (+ 0.6 s.d.) Rod mount: 
Range: <0.1 to 5.9 min Fixed: 73.1% 


Portable: 26.7% 


Trip Length Hull Construction Number of reels/set: Mainline material: Sea State: 
Avg: 4.8 days (+ 3.6 s.d.) Fiberglass: 89% Avg: 3.3 (+ 1.4 s.d.) Monofilament (76.8%), Cable 0 to 2 foot seas: 59% 


Range: 1 to 17 days Wood: 5% Range: 1 to 14 (13.7%), Mono/nylon/poly (3.2%), 3 to 5 foot seas: 31% 
Steel: 4% Other (6.3%) 6 to 8 foot seas: 8% 
Fiberglass/wood: 1% Test: 8+ foot seas: 2% 
Unknown: 1% Avg: 258.3 Ibs (+ 233.6 s.d.) 


Range: 12 to 1,400 Ibs 


’ Crew size Engine Horsepower Hooks: Subline material: Bottom type: 
0 to 4 individuals Avg: 326.9 hp Avg: 26.1 hooks (+ 44.8 s.d.) Monofilament: 97.8% Rock: 67% 
(excluding captain) (+ 195.6 s.d.) Range: 1 to 330 hooks Test: Unknown: 16% 
Range: 40 to1200 hp Type: Circle hooks (83.3%), Avg: 127.2 Ibs (+ 58.5 s.d.) Shell: 2% 
J-hooks (12.7%), double J-hooks Range: 10 to 800 Ibs Coral: 4% 
(3.1%), other (0.8%) Mud: 5% 
Size: 8 aught (43%), 9 aught (20%) Sand: 5% 
Range: 1 to 18 aught Wreck: 1% 
Brand: Mustad® (44%), Eagle 
Claw® (0.4%) Hooks/Reel: Fishing State: 
Avg: 7.4 hooks (+ 10.8 s.d.) On anchor: 68% 
Range: 1 to 45 hooks Drifting: 24% 
Trolling: 2% 


Unknown: 6% 


Marine Fisheries Review 


| 
¢ 
H 
i 
6 


Red Grouper Disposition 
and Size Composition 


All 40,992 red grouper caught using 
longline were in the eastern Gulf of 
Mexico, with the exception of two 
individuals recorded in the western 
Gulf. Based on visual observations, 
the majority (43%) of the fish were 
released alive, 40% were kept, 12% 
were discarded dead, and 6% were of 
unknown condition.'’ One red grouper 
was used for bait. 

A total of 36,764 red grouper were 
measured and ranged from 4 to 37 in 
TL with the mode of 4,440 individuals 
at 18 in TL (Fig. 3). Of these, 32% of 
the fish caught were <18 in TL, the legal 
minimum size, with 69% released alive, 
19% discarded dead, 11% discarded in 
an unknown condition, and 0.3% kept. 
Of the 68% of red grouper >18 in TL, 
62% were kept, 26% were released 
alive, 8% were discarded dead, and 3% 
discarded in an unknown condition. 

Depths of red grouper captures ranged 
from 19.3 to 120.5 fm. Most (67%) red 
grouper were caught between 20—25 fm, 
followed by 26-30 fm (21%), 31-35 fm 
(5%), and 36-40 fm (4%). Catch was 
<1% for the remaining zones (Fig. 4). 


CPUE and Discard CV 


Mean CPUE for all species and 
dispositions combined was 0.0095 fish 
per hook-hour (+ 0.0002 SE; Table 5). 
The catch rate estimate for red grouper 
was 0.0021 fish kept per hook-hour 
(+ 0.0001 SE). Spatial CPUE density 
(numbers of fish kept per 1,000 hook- 
hour) for dominant species for all years 
combined is depicted (Fig. 5-9). Red 
grouper were caught and retained pri- 
marily in statistical areas 2 through 8, 
with highest density CPUE observed in 
Statistical area 5. 

A similar pattern was detected for 
blueline tilefish with highest density 
CPUE in the eastern Gulf of Mexico. 
Yellowedge grouper, tilefish, and scamp, 
Mycteroperca phenax, were distributed 
throughout the Gulf with high CPUE 
observed in deeper waters of the west- 
ern Gulf. Clusters of significantly high 


'7Percentages may not equal 100% due to rounding. 
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Table 4.—Number of fish observed using longline (n=1,503 sets) and vertical line (n=7,391 sets) gear in the Gulf of 
Mexico from July 2006 to December 2009. 


Common name 


Scientific name 


Longline 


Vertical line 


Total 


Red grouper 

Red snapper 

Vermilion snapper 
Yellowedge grouper 

Red porgy 

Blueline tilefish 

Gag 

Tilefish 

Atlantic sharpnose shark 
Scamp 


King snake eel 
Smooth dogfish 
Sharks grouped 
Snowy grouper 
Gray snapper 
King mackerel 
Greater amberjack 
Blacknose shark 
Gray triggerfish 
Chub mackerel 


Yellowtail snapper 
Pinfish 

Blue runner 
Speckled hind 
Lane snapper 
Tomtate 

Almaco jack 
Knobbed porgy 
Spotted hake 
Palespotted eel 


Jolthead porgy 

Mutton snapper 
Sharksucker 

Banded rudderfish 
White grunt 

Little tunny 

Lesser amberjack 
Southern hake 
Spinycheek scorpionfish 
Great barracuda 


Nurse shark 
Sand perch 

Gulf hake 

Silky shark 
Lemon shark 
Bearded brotula 
Dolphin 
Blackedge moray 
Blacktail moray 
Moray (genus) 


Warsaw grouper 
Jack (genus) 
Blacktip shark 
Black sea bass 
Remora 

Florida pompano 
Tiger shark 

Spotted moray 
Creole-fish 
Purplemouth moray 


Black grouper 

Cobia 

Sand seatrout 

Leopard toadfish 
Dogfish (genus) 

Bank sea bass 

Bluefish 

Scalloped hammerhead 
Cubera snapper 
Dogfish 


Whitebone porgy 
Inshore lizardfish 
Crevalle jack 


Epinephelus morio 

Lutjanus campechanus 
Rhomboplites aurorubens 
Epinephelus flavolimbatus 
Pagrus pagrus 

Cauilolatilus microps 
Mycteroperca microlepis 
Lopholatilus chamaeleonticeps 
Rhizoprionodon terraenovae 
Mycteroperca phenax 


Ophichthus rex 
Mustelus canis 

General sharks 
Epinephelus niveatus 
Lutjanus griseus 
Scomberomorus cavalla 
Seriola dumerili 
Carcharhinus acronotus 
Balistes capriscus 
Scomber japonicus 


Ocyurus chrysurus 
Lagodon rhomboides 
Caranx crysos 

Epinephelus drummondhayi 
Lutjanus synagris 
Haemulon aurolineatum 
Seriola rivoliana 

Calamus nodosus 
Urophycis regia 

Ophichthus puncticeps 


Calamus bajonado 
Lutjanus analis 

Echeneis naucrates 
Seriola zonata 

Haemulon plumieri 
Euthynnus alletteratus 
Seriola fasciata 
Urophycis floridana 
Neomerinthe hemingwayi 
Sphyraena barracuda 


Ginglymostoma cirratum 
Diplectrum formosum 
Urophycis cirrata 
Carcharhinus falciformis 
Negaprion brevirostris 

Brotula barbata 

Coryphaena hippurus 
Gymnothorax nigromarginatus 
Gymnothorax kolpos 
Gymnothorax sp 


Epinephelus nigritus 
Seriola sp. 
Carcharhinus limbatus 
Centropristis striata 
Remora remora 
Trachinotus carolinus 
Galeocerdo cuvier 
Gymnothorax moringa 
Paranthias furcifer 
Gymnothorax vicinus 


Mycteroperca bonaci 
Rachycentron canadum 
Cynoscion arenarius 
Opsanus pardus 
Squalus 

Centropristis ocyurus 
Pomatomus saltatrix 
Sphyrna lewini 
Lutjanus cyanopterus 
Mustelus sp. 


Calamus leucosteus 
Synodus foetens 
Caranx hippos 


40,992 
2,456 
139 
6,983 
568 
3,591 
723 
2,199 
2,142 
993 


1,573 
1,284 
1,025 
949 
110 
16 
270 
816 
29 


13,855 
27,669 
26,045 
104 
6,120 
23 
2,624 
45 


54,847 
30,125 
26,184 
7,087 
6,688 
3,614 
3,347 
2,244 
2,225 
1,995 


1,585 
1,319 
1,121 
1,117 
932 
902 
883 
848 
837 
818 


781 
599 
532 
523 
509 
495 
492 
408 
380 
288 


286 
285 
277 
267 
263 
255 
239 
230 
211 
198 


197 
168 
168 
166 
165 
161 
158 
149 
147 
141 


134 
132 
127 
127 
117 
116 
113 
112 
107 
106 


101 


continued 


1,002 
35 
96 
168 
822 
886 
613 
32 4 
808 
0 818 as 
11 770 
1 598 
7 525 
492 31 
93 416 3 
1 494 
39 453 
12 396 iy 
377 3 
288 0 te 
132 154 
265 20 : 
213 64 pe 
12 255 “ae 
4 259 
127 128 
20 219 
230 0 
208 3 
153 45 a 
163 34 
38 130 
168 
95 71 4 
157 8 
148 13 
91 67 
141 8 
144 3 i 
133 8 4 
80 54 
114 18 
87 40 : 
0 127 
37 80 
2 114 
107 6 
83 29 4 
0 107 
97 9 
67 34 = 
1 72 28 100 2 
24 74 98 : 
4 79 13 92 “yf 
92 0 92 5 
20 61 81 
2 78 80 
76 2 78 ; 
76 2 78 
i 72 5 77 te 
6 67 73 ; 
9 59 68 
} 
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c Table 4.—(Continued). CPUE for red grouper were located in 
Common name Scientific name Longline Vertical line Total statistical areas 3 through 8 (Fig. 10). 
: x pe For all kept species, clusters of sig- 
Queen snapper Etelis oculatus 1 
Red drum 22 43 65 nificantly high CPUE were detected in 
Grunt (genus) — = statistical areas 5, 14, 15, and 16 (Fig. 
: Sandbar shark Carcharhinus plumbeus 59 2 61 11): Highest discard CPUE was evident 
lizardfish in statistical areas 3 through 6 (Fig. 12). 
jar jack aranx ruber 
CV estimates (Table 7) for discarded 
Blackfin tuna Thunnus atlanticus 49 9 58 ar 
3 Blackbelly rosefish Helicolenus dactylopterus 42 10 52 red grouper, red — greater amber 
Cuban dogfish Squalus cubensis 49 1 50 jack, Seroila dumerili; and gag, Mycte- 
Clearnose skate Raja eglanteria 50 0 50 r(< 
Wenchman Pristipomoides aquilonaris 23 25 4g Toperca mic rolepis, were low (s0.1 
: Smalltail shark Carcharhinus porosus 48 0 48 Several other species of grouper; jacks, 
44 king mackerel, Scomberomorus cavalla; 
: Bull shark Carcharhinus leucas 43 0 43. and cobia, Rachycentron canadum, had 
y Silver seatrout Cynoscion nothus 20 18 38 values <0.5 
Lizardfish (family) Synodontidae 31 5 36 : 
5 Gulper shark Centrophorus granulosus 35 0 35 Vertical Line 
Sharpnose sevengill shark Heptranchias perlo 33 0 33 
Spinner shark Carcharhinus brevipinna 28 2 30 Allocation of Sampling Effort 
Sand diver Synodus intermedius 27 2 29 ‘ ies 
Bigeye 0 29 29 Data from 237 trips were collected 
aboard 157 vertical line vessels from 
ittlehead porgy alamus prorigens 
Gulf toadfish Opsanus beta 21 4 25 July 2006 through December 2009, with 
Great hammerhead Sphyrna mokarran 2A 0 24 atotal of 89.015 fish processed (Tables 3 
Chain dogfish Scyliorhinus retifer 24 0 24 and 4). Locations for 7,384 vertical line 
Short bigeye Pristigenys alta 3 20 23 oe eles io 
Ocean triggerfish Canthidermis sufflamen 0 23 23 sets are depicted (Fig. 13). Effort data 
Squirrelfish Holocentrus adscensionis 3 19 22 (5.266 h; 190.202 hooks) were available 
Cubbyu Pareques umbrosus 0 22 22 
99 for 7,285 sets. Approximately 37% of 
A Night shark Carcharhinus signatus 20 0 20 the sampled reels had no catch reported 
Yellowmouth grouper Mycteroperca interstitialis 9 10 19 
‘ Triggerfish (family) Balistidae 0 19 19 during a set. The majority (75% ) of sets 
Rock hind Epinephelus adscensionis 1 18 19 were in the eastern Gulf of Mexico. 
7 12 19 However, the highest concentrated effort 
Wahoo Acanthocybium solandri 10 8 18 (74%), based on hook-hours, occurred 
Reticulate moray Muraena retifera 18 0 18 
Siarkher drum Equetus iwamotoi 0 18 ig inthe western Gulf of Mexico (Fig. 14). 
Round scad Decapterus punctatus = ue - By season, 23% of the effort occurred 
Hake (genus) Urophycis sp is 
Jack (family) Carangidae 4 12 16 from January through March; 21% 
Graysby onsentete April through June; 33% July through 
Tattler Serranus phoebe 0 14 14 jee 096 = 
Saquirrelfishes (family) Holocentridae 3 "1 14 September; and 22% October through 
~ 7 ; 
Rainbow runner Elagatis bipinnulata 6 8 14 December for all years combined. 
Black margate Anisotremus surinamensis 14 0 14 
Bigeye scad Selar crumenophthalmus 0 14 14 Species Composition 
Biuntnose sixgill shark Hexanchus griseus 13 0 13 : : 
Red hind Epinephelus guttatus 2 11 13 Of the 89.015 fish (178 taxa) sampled, 
71% of the individuals were kept, 19% 
corpionfish Scorpaena sp 9 3 12 i i 
Rock sea bass Centropristis philadelphica 8 4 12. were released alive, 6% were discarded 
: Horse-eye jack Caranx latus 0 12 12 Oo 
: Toadfish (genus) Opsanus sp 1 0 +; dead : 1 7o were discarded in an unknown 
: condition, and 4% were retained for 
; Silk snapper Lutjanus vivanus 7 4 : x : 
Longiall bess nuaniidietnne 10 bait (Tables 5 and 8). By number, red 
Dusky shark Carcharhinus obscurus "1 0 snapper ranked highest in catch com- 
Bigeye sixgill shark Hexanchus nakamurai 11 0 $42 +1: 
Atlantic croaker Micropogonias undulatus 0 11 position at 3 1%. Vermilion snapper 
Smooth putter Lagocephalus laevigatus 10 0 comprised 29% of the catch, followed 
Largescale lizardfish Saurida brasiliensis 9 0 b d 16% d P 5 
Atlantic spadefish Chaetodipterus faber 0 9 y re grouper ( o ), re porgy, agrus 
Hardhead catfish Arius felis 0 8 pagrus (7%); gag (3%), and the remain- 
Grunt (family) Haemulidae 8 0 ‘no species combined ( 14% ) 
0). 
Goldface tilefish Caulolatilus chrysops 1 7 8 Vermilion snapper, red snapper, red 
Southern stingray Dasyatis americana 6 1 7 . 
Cusk-eel (family) Ophidiidae 5 2 7 grouper, and red porgy comprised 86% 
Barracuda (genus) Sphyraena sp 6 1 7 H Wi H 
Atlantic cutlassfish Trichiurus lepturus 2 5 7 of the 63,351 individuals in the kept 
Spiny dogfish Squalus acanthias 6 0 6 category. Three species (red snapper, 
°¢ red grouper, and vermilion snapper) 
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accounted for 80% of the released alive 
category. Of the 16,872 individuals re- 
leased alive, 35% exhibited visual signs 
of stress, while 61% exhibited a normal 
appearance. 

Of the 2,805 individuals used for 
bait, the species caught and used most 
often were chub mackerel, Scomber 
japonicus (29%); pinfish, Lagodon 
rhomboides (20%); and tomtate, Hae- 
mulon aurolineatum (16%). Red snap- 
per, vermilion snapper, and red grouper 
comprised 87% of 5,185 individuals in 
the discarded dead category. Minimum 
assumed mortality for these species was 
approximately: red snapper (28%), ver- 
milion snapper (41%), and red grouper 
(11%). The fate of 802 individuals was 
not determined. 


Red Snapper Disposition 
and Size Composition 


A total of 27,669 red snapper were 
sampled on vertical line gear. Statisti- 
cal areas of capture ranged from 3 to 
21, with no reported takes in statistical 
area 12. Approximately 77% of the red 
snapper were captured in the western 
Gulf of Mexico, with the remaining 23% 
captured in the eastern Gulf. The major- 
ity (65%) of the fish were kept. Based on 
visual observations, 24% were released 
alive, 10% were discarded dead, and 
1% discarded in an unknown condition. 

A total of 25,650 red snapper were 
measured and ranged from 6 to 41 in 
TL, with the mode of 4,102 individuals 
at 15 in TL(Fig. 15). Of these, 92% were 
>13 in TL, the legal minimum size. Ap- 
proximately 8% were <13 in TL, with 
31% of the individuals discarded dead. 

Depths of red snapper capture ranged 
from 3.3 to 305 fm. Most (29%) red 
snapper were caught in waters less than 
20 fm, followed by 20-25 fm (26%), 
and 31-35 and 26-30 fm (13% each; 
Fig. 16). The remaining depth zones 
comprised 19%. No depth values were 
recorded for 762 red snapper. 


CPUE and Discard CV 


Mean CPUE for all species and dis- 
positions was 0.9369 fish per hook-hour 
(+ 0.0311 SE; Table 5). Red snapper 
mean catch rate was 0.2214 fish kept 
per hook-hour (+ 0.0150 SE). Spatial 


Table 4.—(Continued). 


Common name 


Scientific name 


Longline 


Vertical line 


Total 


Shortfin mako 
Margate 
Grass porgy 
Atlantic bonito 


Swordfish 
Sailors choice 
Honeycomb moray 


Hammerhead (genus) shark 


Green moray 

Florida smoothhound 
Finetooth shark 
Thresher shark 
Atlantic stingray 
Starfish (family) 


Spider (genus) crab 
Southern flounder 
Snake eel (family) 
Sea bass (family) 
Sailfish 

Queen triggerfish 
Puffer (family) 
Porgy (genus) 
Pigfish 

Black snapper 


Anchor tilefish 
Spottail pinfish 
Spanish flag 

Shoal flounder 
Saucereye porgy 
Octopus (genus) 
Guaguanche 
Conger eel (family) 
Conger eel 
Bonnethead 


Black jack 

Black drum 
Bermuda chub 
Yellowfin grouper 
Yellow conger 
Spotfin hogfish 
Southern puffer 
Smooth butterfly ray 
Pufferfish (genus) 
Porgie (family) 


Oyster toadfish 
Mackerel (family) 
Lefteye flounder (family) 
Fish (superclass) 
Dusky flounder 

Drum (family) 

Cero 

Broad flounder 

Atlantic angel shark 
Yellow jack 


Whitespotted soapfish 
Threadtail conger 
Stingray (genus) 
Stingray (family) 
Spotted snake eel 
Spanish sardine 
Spanish hogfish 
Skipjack tuna 

Skate (genus) 

Shrimp eel 


Sand tiger 

Saddled grenadier 
Roughtongue bass 
Rosette skate 
Porkfish 

Offshore hake 
Octopus (order) 
Ocellated frogfish 
Marbled grouper 
Mantis (genus) shrimp 


Isurus oxyrinchus 
Haemulon album 
Calamus arctifrons 
Sarda sarda 


Xiphias gladius 
Haemulon parra 
Gymnothorax saxicola 
Sphyrna sp 
Gymnothorax funebris 
Mustelus norrisi 
Carcharhinus isodon 
Alopias vulpinus 
Dasyatis sabina 
Astropectinidae 


Libinia sp. 

Paralichthys lethostigma 
Ophichthidae 
Serranidae 

Istiophorus platypterus 
Balistes vetula 
Tetraodontidae 

Calamus 

Orthopristis chrysoptera 
Apsilus dentatus 


Caulolatilus intermedius 
Diplodus holbrooki 
Gonioplectrus hispanus 
Syacium gunteri 
Calamus calamus 
Octopus sp. 

Sphyraena guachancho 
Congridae 

Conger oceanicus 
Sphyrna tiburo 


Caranx lugubris 
Pogonias cromis 
Kyphosus sectatrix 
Mycteroperca venenosa 
Hildebrandia flava 
Bodianus pulchellus 
Sphoeroides nephelus 
Gymnura micrura 
Sphoeroides sp 
Sparidae 


Opsanus tau 

Scombridae 

Bothidae 

Pisces 

Syacium papillosum 
Sciaenidae 
Scomberomorus regalis 
Paralichthys squamilentus 
Squatina dumeril 

Caranx bartholomaei 


Rypticus maculatus 
Uroconger syringinus 
Dasyatis sp 
Dasyatidae 
Ophichthus ophis 
Sardinella aurita 
Bodianus rufus 
Katsuwonus pelamis 
Raja 

Ophichthus gomesi 


Carcharias taurus 
Caelorinchus caelorhincus 
Holanthias martinicensis 
Raja garmani 
Anisotremus virginicus 
Merluccius albidus 
Octopoda 

Antennarius ocellatus 
Epinephelus inermis 
Squilla sp. 
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Longline 2006 - 2009 


0.591 - 155.563 


GNM 155.564 - 310.535 
MM 310.536 - 465.508 


Gulf of Mexico 


Sum of Effort - (Hooks*Hours)*1000 
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Figure 2.—Distribution of sampling effort (hook-hours) based on observer coverage of the U.S. Gulf of Mexico bottom 
longline reef fish fishery from August 2006 through November 2009. 


CPUE density (numbers of fish kept 
per 1,000 hook-hours) for dominant 
species caught using vertical line gear 
is depicted in Figures 17 through 21. 
Red snapper were caught and retained 
throughout the Gulf, with highest 
density CPUE observed in statistical 


Table 4.—(Continued). 


area 11. Similarly, vermilion snapper 
occurred in both Gulf regions with a 
spatial density similar to red snapper. 
Red grouper were concentrated in the 
eastern Gulf, with the highest CPUE 
density observed in statistical areas 3,4, 
and 8. High density CPUE for red porgy 


Common name 


Scientific name 


Vertical line 


Lookdown 

Longspine squirrelfish 
Jack (genus) 

Gulf hagfish 

Gulf flounder 
Gafftopsail catfish 
Dog snapper 
Decapod (order) 

Big roughy 

Cusk-eel (genus) 


Cownose ray 
Cottonwick 
Cottonmouth jack 
Cardinal soldierfish 
Butterfly ray 
Biuntnose stingray 
Blackline tilefish 
Bigeye tuna 
Barrelfish 

Bank cusk-eel 
Atlantic moonfish 


Total 


Selene vomer 
Holocentrus rufus 
Caranx 

Eptatretus springeri 
Paralichthys albigutta 
Bagre marinus 
Lutjanus jocu 
Decapoda 
Gephyroberyx darwinii 
Lepophidium 


Rhinoptera bonasus 
Haemulon melanurum 
Uraspis secunda 
Plectrypops retrospinus 
Gymnura sp 

Dasyatis say 
Cauloiatilus cyanops 
Thunnus obesus 


Hyperoglyphe perciformis 


Ophidion holbrooki 
Selene setapinnis 


was found primarily in the eastern Gulf, 
with the exception of statistical area 16. 
Gag were caught and retained primarily 
off Florida, predominantly in statistical 
areas 5-8. 

Cluster locations of statistically 
significant high CPUE for retained 
red snapper were most pronounced in 
statistical areas 8 through 14, 16, and 
17 (Fig. 22). For all retained species, 
clusters of significantly high CPUE 
were detected primarily in the western 
Gulf (Fig. 23). Conversely, highest 
discard CPUE values were observed 
in the eastern Gulf in statistical areas 5 
through 7 (Fig. 24). 

Based on number discarded, CV es- 
timates for Federally managed species 
caught on vertical line gear (Table 9) 
were low for red grouper, red snapper, 
vermilion snapper, gag, and greater am- 
berjack (<0.1). Several other species of 
grouper, jacks, gray triggerfish, Balistes 
capriscus; king mackerel, and red drum, 
Sciaenops ocellatus, had values less than 
or equal to 0.5. Higher CV estimates for 
other species of importance, including 
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several species of snapper and grouper, 
were detected. 


Interactions with 
Protected Species in 
the Reef Fish Fishery 


Twenty sea turtles were captured on 
observed trips utilizing longline gear 
from 2006 to 2009; three occurred 
during the electronic monitoring pilot 
project. One sea turtle was captured on 
vertical line gear (bandit) during the 
same time period. Sea turtle mortality 
and projected take estimates by gear 
type were reported by SEFSC.!$ 


Discussion 


To gain a greater understanding 
of catch rates, bycatch composition, 


'8SEFSC. 2009. Estimated takes of sea turtles 
in the bottom longline portion of the Gulf of 
Mexico reef fish fishery July 2006 through 
December 2008 based on observer data. U.S. 
Dep. Commer., NOAA, NMFS Southeast Fish. 
Sci. Cent. Contrib. PRD-08/09-07, March 
2009, 23 p. [Updated 4/2009, Erratum: updated 
6/2009]. 
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Table 5.—Species composition and disposition by gear type observed from July 2006 to December 2009. 
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Figure 3.—Size and fate of red grouper caught on bottom longline gear based on observer coverage of the U.S. Gulf 
of Mexico reef fish fishery from August 2006 through November 2009. 
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Vertical line 


73,205 fish of 183 taxa 


89,015 fish of 178 taxa 


Kept: 46% 
Red grouper: 49% 
Yellowedge grouper: 21% 
Tilefish: 6% 
Blueline tilefish: 5% 


Released alive: 35% 
(42% stressed: air bladder expansion and/or eyes 
protruding; 46% normal; 12% not recorded) 
Red grouper: 69% 
Atlantic sharpnose shark, Smooth dogfish, Red snapper: 
5% each 


Discarded dead: 12% 
Red grouper: 54% 
Blueline tilefish: 15% 
Atlantic sharpnose shark: 8% 
Red snapper: 5% 


Unknown: 4% 
Red grouper: 77% 
Atlantic sharpnose shark, Gulf hake, Grouped sharks: 
3% each 


Kept for bait: 3% 
King snake eel: 29% 
Palespotted eel: 11% 
Little tunny: 5% 


Mean CPUE (fish/hook hour): 
All: 0.0095 (+ 0.0002) 
Kept: 0.0043 (+ 0.0001) 
Red grouper: 0.0021 (+ 0.0001) 


Sea turtle captures: 19 


Kept: 71% 
Vermilion snapper: 37% 
Red snapper: 28% 
Red grouper: 12% 
Red porgy: 9% 


Released alive: 19% 


(35% stressed: air bladder expansion and/or eyes 
protruding; 61% normal; 4% not recorded) 


Red snapper: 39% 
Red grouper: 34% 
Vermilion snapper: 7% 


Discarded dead: 6% 
Red snapper: 53% 
Vermilion snapper: 21% 
Red grouper: 13% 


Unknown: 1% 
Vermilion snapper: 45% 
Red snapper: 43% 
Red grouper: 5% 


Kept for bait: 4% 
Chub mackerel: 29% 
Pinfish: 20% 
Tomtate: 16% 


Mean CPUE (fish/hook hour): 
All: 0.9369 (+ 0.0311) 
Kept: 0.6500 (+ 0.0221) 
Red snapper: 0.2214 (+ 0.0150) 


Sea turtle captures: 1 


aby 
4500 
| | 
4000 : 
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Table 6.—Number, condition (when brought onboard), and fate of fish Species with n>25 caught using longline gear in the Gulf of Mexico from August 2006 to November 2009. 


12 


Fate upon release Kept Released alive Kept for bait Discarded dead Unknown 
Condition upon capture Live Live Live Live Live 
> > > £ 
Common name 
Red grouper 40,992 16413 4186 10,402 259 17,475 5,078 9,543 1 1 4843 1,010 2811 760 2,260 98 890 
Yellowedge grouper 6.983 6,932 251 5,759 918 5 4 6 1 5 15 4 11 25 4 12 
Blueline tilefish 3,591 1,767 551 1,179 37 417 152 264 67 43 14 10 1,331 212 782 332 9 3 5 
Red snapper 2,456 784 501 269 3 1,161 376 702 1 1 450 132 208 92 60 16 35 
Tilefish 2,199 2.130 996 1,036 93 9 8 1 J 3 1 32 6 10 16 24 3 21 
Atlantic sharpnose shark 2,142 20 12 1 7 1,280 1,264 4 50 35 1 14 699 280 2 379 93 79 
King snake ee! 1,573 2 2 714 711 1 692 672 4 5 150 110 11 8 15 6 
Smooth dogfish 1,284 1 1 1,176 1,173 2 52 52 44 31 10 11 8 1 
Sharks grouped 1,025 1 1 710 701 13 13 275 141 129 26 10 
Scamp 993 955 453 439 14 22 10 5 13 3 6 a 3 1 
Snowy grouper 949 941 114 771 55 2 1 1 6 1 2 2 
Blacknose shark 816 6 6 576 572 15 9 6 162 92 58 57 54 
Gag 723 673 234 417 41 14 22 7 1 4 2 2 1 
Red porgy 568 507 363 119 2 16 13 3 29 24 2 1 10 3 4 3 6 6 
Speckled hind 492 453 99 324 28 17 5 9 22 7 4 
Spotted hake 377 7 3 4 2 2 68 2 60 6 262 163 99 38 5 32 
Palespotted ee! 288 9 7 271 261 1 6 d 1 2 1 
Greater amberjack 270 124 112 1 7 99 97 14 14 22 13 1 8 11 8 
Mutton snapper 265 264 216 47 1 1 1 
Southern hake 230 7 2 5 5 3 2 50 6 37 6 135 a 116 15 33 2 31 
Sharksucker 213 1 1 148 128 47 47 5 4 12 1 
Spinycheek scorpionfish 208 202 62 114 25 5 1 3 1 1 1 
Gulf hake 168 13 4 8 2 2 65 56 9 88 84 
Nurse shark 163 142 127 1 20 11 
Lemon shark 157 153 153 4 1 3 
Great barracuda 153 11 11 15 14 107 79 13 14 7 7 6 5 
Bearded brotula 148 128 81 35 12 1 1 2 1 1 16 1 15 1 1 
Biacktail moray 144 11 11 89 85 4 44 42 2 
Biackedge moray 141 1 1 37 37 81 66 15 16 10 5 6 3 
Vermilion snapper 139 84 18 33 4 32 22 1 11 6 4 11 2 3 4 1 
Moray (genus) 133 9 9 100 78 21 18 5 9 6 1 
Jolthead porgy 132 127 115 3 1 1 1 4 4 
Little tunny 127 1 1 113 14 93 13 2 10 
Jack (genus) 114 71 69 1 5 5 38 = 38 
Gray snapper 110 105 25 49 1 3 2 
Tiger shark 107 97 94 1 1 a 1 1 5 2 
Purplemouth moray 97 4 4 64 47 17 29 15 12 
Silky shark 95 58 57 2 1 1 34 9 24 1 1 
Lane snapper 93 75 18 49 3 7 3 2 1 5 1 2 2 5 
Dogfish (genus) 92 52 52 38 38 2 2 
Dolphin 91 89 22 67 1 1 1 
Blacktip shark 87 7 4 3 55 54 7 5 2 17 1 15 1 1 
Spotted moray 83 19 19 54 27 23 10 3 Fd 
Warsaw grouper 80 78 6 71 1 1 1 1 1 
Leopard toadfish 79 35 20 14 34 18 16 8 5 3 2 1 1 
Cubera snapper 76 76 75 1 
Scalloped hammerhead 76 1 1 56 54 1 1 13 13 5 2 
Dogfish 72 69 68 1 1 1 e 2 
Cobia 72 38 34 1 29 28 2 2 3 3 
Black grouper 67 65 31 15 2 1 
Inshore lizardfish 66 20 3 40 32 1 a 5 1 1 1 
Sandbar shark 59 57 54 2 2 
Clearnose skate 50 9 7 41 39 2 
Cuban dogfish 49 36 36 8 8 5 5 
Blackfin tuna 49 38 17 21 2 2 6 6 2 1 1 1 
Smalitail shark 48 48 48 
Snakefish 44 8 2 33 21 1 11 3 , 
Bull shark 43 42 42 1 1 
Blackbelly rosefish 42 12 11 1 12 9 3 18 2 16 
Offshore lizardfish 41 7 7 26 11 1 13 8 3 3 
Almaco jack 39 19 19 3 3 11 11 6 6 
Sand perch 38 12 5 1 24 18 2 2 1 1 1 
Remora 37 34 34 3 2 
Gulper shark 35 30 30 5 5 
Sevengill shark 33 25 25 8 8 
Lizardfish (family) 31 5 5 23 12 11 2 2 1 
Gray triggerfish 29 26 16 8 3 1 
Spinner shark 28 2 2 15 15 9 8 1 2 
Sand diver 27 25 22 3 2 2 
Total (all species) 73,205 33,335 8,778 21,183 1,583 25,471 11,744 10,628 2.414 1,849 178 320 9,037 2.235 4,258 2,149 2,948 407 1,132 
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Figure 4.—Number of red grouper by size and depth zone caught on bottom longline gear based on observer 
coverage of the U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
92° W 88° W 84° W 
Longline 2006 - 2009 . 
Red Grouper (Epinephelus morio) eres x 
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| 0.030 - 3.348 
HME 3.349 - 6.667 
ces - 9.986 
Figure 5.—CPUE density surface for red grouper kept in the bottom longline fishery based on observer coverage 
of the U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 6.—CPUE density surface for blueline tilefish kept in the bottom longline fishery based on observer coverage of 
the U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 7.—CPUE density surface for yellowedge grouper kept in the bottom longline fishery based on observer coverage 
of the U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 8.—CPUE density surface for tilefish kept in the bottom longline fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 9.—CPUE density surface for scamp kept in the bottom longline fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 10.—Hot Spot Analysis for all kept red grouper in the bottom longline fishery based on observer coverage of the 
U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 11.—Hot Spot Analysis for all kept species in the bottom longline fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 12.—Hot Spot Analysis for all discarded species in the bottom longline fishery based on observer coverage of the 
U.S. Gulf of Mexico reef fish fishery from August 2006 through November 2009. 
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Figure 13.—Distribution of sampling effort (sets) based on observer coverage of the U.S. Gulf of Mexico vertical line 
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and discard mortality associated with 
commercial fishing operations in the 
U.S. Gulf of Mexico reef fish fishery, 
a mandatory observer program was 
established in 2006 based on a pro- 
portional randomized sampling design 


stratified by season, gear, and region. 
Historically, these data, critical for 
population assessments, have not been 
available due to lack of time series and 
limited geographic ranges for affected 
species. 


Table 7.—Coefficient of variation (CV) for Federally managed discarded species caught aboard longline vessels in 
the Gulf of Mexico from August 2006 to November 2009. 


Common name Scientific name n CV 


Red grouper Epinephelus morio 24,081 <0.1 
Red snapper Lutjanus campechanus 1,657 0.1 
Blueline tilefish Caulolatilus microps 1,824 0.1 
Greater amberjack Seriola dumerili 133 0.1 
Gag Mycteroperca microlepis 48 0.1 
Vermilion snapper Rhomboplites aurorubens 43 0.2 
Tilefish Lopholatilus chamaeleonticeps 67 0.2 
Cobia Rachycentron canadum 27 0.2 
Speckled hind Epinephelus drummondhayi 39 0.2 
Yellowe ol Epinephelus flavolimbatus 50 0.2 
Lesser amberjack 19 0.3 
Lane snapper 18 0.3 
es aquilonaris 17 0.3 

Epinephelus niveatus 8 0.4 

Mycteroperca phenax 37 0.4 

Scomberomorus cavalla 6 0.4 

5 0.5 

10 0.5 

16 0.6 

2 0.7 

2 0.7 

Balistes capriscus 2 0.7 

Mycteroperca bonaci 2 0.7 

s Ocyurus chrysurus 3 0.7 

Mutton snapper Lutjanus analis 1 1.0 
Rock hind Epinephelus adscensionis 1 1.0 


Data from this observer program re- 
vealed relatively high species richness 
from the two primary gears (longline n= 
183 taxa; and vertical line n = 178 taxa). 
While diversity was high, red grouper 
and yellowedge grouper (in longline), 
and red snapper and vermillion snapper 
(in vertical line), comprised more than 
60% by number of the species caught. 
These findings are similar to those de- 
scribed by Stephen and Harris (2010) of 
the snapper-grouper vertical line fishery 
off South Carolina. Their data revealed 
high overall diversity; however, a small 
number of species (17) accounted for 
90% of catch. 

Hale et al. (2010), through a man- 
datory bottom longline observer pro- 
gram, examined species composition 
and disposition of fish captured from 
longline sets targeting reef fish in the 
Gulf of Mexico and found, in order of 
abundance, that red grouper, blueline 
tilefish, tilefish, and yellowedge grouper 
comprised 76% of catch. In our current 
study, these four species accounted for 
73% of the catch captured on longline 
gear. Moreover, disposition of these 
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cal line reef fish fishery from July 2006 through December 2009. 


Figure 14.— Distribution of sampling effort (hook-hours) based on observer coverage of the U.S. Gulf of Mexico verti- 
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species was similar between the two 
programs for red and yellowedge grou- 
per. Blueline tilefish and tilefish discard 
proportion rates were more variable, and 
most likely related to the 15 May 2009 


tilefish quota closure. 


In our current study, 46% of the 
individuals, predominately red and 
yellowedge grouper, were kept in 
longline. In vertical line, a larger per- 
centage (71%) was kept and comprised 
primarily of vermilion and red snapper. 


While species-specific minimum size 
limits differ by region, Rudershausen 
et al. (2007), Stephen and Harris 
(2010), and Scott-Denton® reported 
low discard proportions for the verti- 
cal line trips; however, low discard 


Table 8.--Number, condition (when brought onboard), and fate of fish species with n>25 caught using vertical line gear in the Gulf of Mexico from July 2006 to December 2009. 


Fate upon release Kept 


Released alive 


Kept for bait 


Discarded dead Unknown 


Condition upon capture 


Common name 


Live 


Live 


Live Live 


@ 


> 


Red snapper 
Vermilion snapper 
Red grouper 

Red porgy 

Gag 

Scamp 

King mackerel 
Gray snapper 
Chub mackerel 
Gray triggerfish 
Yellowtail snapper 
Greater amberjack 
Pinfish 

Biue runner 
Tomtate 

Almaco jack 

Lane snapper 
Knobbed porgy 
White grunt 
Banded rudderfish 
Lesser amberjack 
Snowy grouper 
Joithead porgy 
Sand perch 

Little tunny 

Black seabass 
Florida pompano 
Creole-Fish 
Yellowedge grouper 
Sharks grouped 
Atlantic sharpnose shark 
Remora 

Bluefish 78 
Sand seatrout 74 
Silky shark 71 
Whitebone porgy 67 
Dolphin 67 
Sharksucker 64 
Grunt (genus) 63 
Spanish mackerel! 62 
Bank seabass 61 
Crevalle jack 59 
Bar jack 57 
Warsaw grouper 

Queen snapper 

Sheepshead 

Tilefish 

Great barracuda 

Red drum 

Blacktip shark 

Smooth dogfish 

Nurse shark 

Black grouper 

Blacknose shark 

Speckled hind 

Spotted moray 

Bigeye 

Cobia 

Seatrout (genus) 

Wenchman 


27,669 
26,045 
13,855 
6,120 
2,624 
1,002 
886 
822 
818 
808 
770 722 
613 171 
598 8 
525 129 
494 2 
453 285 
416 388 
396 377 
259 118 
255 
219 
168 
154 
130 
128 
127 
114 
107 
104 
96 
83 
80 


11,368 
21,994 
1,920 
5,022 
874 


23,240 
7,445 
5,971 
1,565 

898 
868 
775 


751 


Total (all species) 89,015 63,351 46,602 
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Figure 15.—Size and fate of red snapper caught on vertical line gear based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 16.—Number of red snapper by size and depth zone caught on vertical line gear based on observer 
coverage of the U.S. Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 17.—CPUE density surface for red snapper kept in the vertical line fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 18.—CPUE density surface for vermilion snapper kept in the vertical line fishery based on observer coverage of 
the U.S. Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 19.—CPUE density surface for red grouper kept in the vertical line fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 20.—CPUE density surface for red porgy kept in the vertical line fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 21.—CPUE density surface for gag kept in the vertical line fishery based on observer coverage of the U.S. Gulf of 
Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 22.—Hot Spot Analysis for all kept red snapper in the vertical line fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 23.—Hot Spot Analysis for all kept species in the vertical line fishery based on observer coverage of the U.S. Gulf 
of Mexico reef fish fishery from July 2006 through December 2009. 
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Figure 24.—Hot Spot Analysis for all discarded species in the vertical line fishery based on observer coverage of the U.S. 
Gulf of Mexico reef fish fishery from July 2006 through December 2009. 
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proportions may still adversely affect 
long-lived stocks. 

Discard mortality rates are highly 
variable and influenced by a number of 
factors, including species-specific life 
history characteristics (Coleman et al., 
2000; Patterson et al., 2002; Nieland et 
al., 2007), season (Render and Wilson, 
1994) depth, and method of capture 
and release (Gitschlag and Renaud, 
1994; Collins et al., 1999, Dorf, 2003: 
Rummer, 2007; Burns et al.’). Using the 
Marine Recreational Fishery Statistic 
Survey data from 1981-99 and findings 
from 53 release mortality studies, Bar- 
tholomew and Bohnsack (2005) found 
significant mortality factors related to 
hook location, bait removal, hook type, 
capture depth, water temperature, and 
handling time. 

Through a tagging study conducted 
off the coast of Alabama, Patterson et al. 
(2002) indirectly estimated discard mor- 
tality of 13.5% for red snapper and <1% 
for gray triggerfish, based on surface re- 
lease observations and recapture rates of 
fish caught with recreational gear. Red 
snapper (<18 in TL) comprised 93% of 
the released fish from a Texas headboat 
survey, of these 60.6% were released 
alive, 22.8% swam erratically, 15.2% 
floated, and 1.4% were discarded dead 
(Dorf, 2003). Diamond and Campbell 
(2009) examined red snapper caught 
on hook and line at three petroleum 
production platforms off south Texas 
and found immediate mortality at 17%; 
however, through the use of an injury 
status condition index, delayed mortality 
was estimated to be 64%. 

Variable minimum assumed mortal- 
ity rates and discard proportions may 
also be attributed to regulatory changes 
in minimum size limits and through 
implementation of IFQ requirements 
for several species, notably, red snap- 
per, red grouper, and tilefish. Minimum 
assumed mortality (all discarded spe- 
cies combined) in this study was 24% 
in longline and 23% in vertical line. By 
species, immediate mortality for red 
grouper was 20% in longline and 11% 
in vertical line, with minimum assumed 
mortality for red snapper of 27% and 
28%, in longline and in vertical line, 
respectively. 
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Table 9.—Coefficient of variation (CV) for Federally-managed discarded species caught aboard vertical line vessels 
in the Gulf of Mexico from July 2006 to December 2009. 


Common name 


Scientific name n CV 
Red grouper Epinephelus morio 6,597 <0.1 
Red snapper Lutjanus campechanus 19,227 <0.1 
Vermilion snapper Rhomboplites aurorubens 5,754 <0.1 
Gag Mycteroperca microlepis 1,096 <0.1 
Greater amberjack Seriola dumerili 621 <0.1 
Lesser amberjack Seriola fasciata 136 0.2 
Gray triggerfish Balistes capriscus 124 0.3 
Warsaw grouper Epinephelus nigritus 32 0.3 
Snowy grouper Epinephelus niveatus 32 0.3 
King mackerel Scomberomorus cavalla 20 0.3 
Banded rudderfish Seriola zonata 363 0.3 
Scamp Mycteroperca phenax 189 0.3 
Cobia Rachycentron canadum 24 0.3 
Goliath grouper Epinephelus itajara 12 0.4 
Speckled hind Epinephelus drummondhayi 24 04 
Yellowedge grouper Epinephelus fiavolimbatus 28 0.4 
Red drum Sciaenops ocellatus 114 0.4 
Lane snapper Lutjanus synagris 79 0.4 
Wenchman Pristipomoides aquilonaris 52 0.4 
Blueline tilefish Caulolatilus microps 8 0.5 
Red hind Epinephelus guttatus 11 0.5 
Rock hind Epinephelus adscensionis 4 0.5 
Yellowtail snapper Ocyurus chrysurus 48 0.6 
Gray snapper Lutjanus griseus 49 0.6 
Spanish mackerel Scomberomorus maculatus 18 0.7 
Black grouper Mycteroperca bonaci 2 0.7 
Queen snapper Etelis oculatus 3 0.7 
Silk snapper Lutjanus vivanus 3 1.0 
Tilefish Lopholatilus chamaeleonticeps 1 1.0 
Mutton snapper Lutjanus analis 1 1.0 
Yellowmouth grouper Mycteroperca interstitialis 1 1.0 


Stephen and Harris (2010) reported 
immediate mortality range of 33—100% 
for vertical line trips targeting vermilion 
snapper off South Carolina, with >90% 
mortality observed for gray triggerfish, 
greater amberjack, scamp, and red snap- 
per. Nieland et al. (2007), using four 
release condition categories, similar but 
more detailed than that of this study, 
assessed the fate of red snapper regula- 
tory discards aboard commercial verti- 
cal line vessels operating primarily off 
Louisiana and found 69% of discarded 
red snapper were either dying or dead 
when released. 

Rudershausen et al. (2007) examined 
discard composition in the commercial 
snapper-grouper fishery in North Caro- 
lina and found low (<10%) immediate 
release morality for vermilion snapper, 
gag, and red grouper; moderate (14%) 
mortality for red porgy; and high (23%) 
immediate mortality for scamp. 

In our study, red snapper ranged from 
6-41 in TL with a mode of 15 in TL. 
Nieland et al. (2007), using specimens 
collected from commercial red snapper 
landings, described a similar unimodal 
distribution with the mode at 400 mm 
(15.7 in) TL, noting that 98% were less 
than 600 mm (23.6 in) TL. Red grou- 


per length frequency data from NMFS 
bottom longline surveys in the Gulf of 
Mexico from 2000 through 2005 de- 
picted a distribution range of approxi- 
mately 10—34 in TL with a mode 18 in 
TL (Ingram et al.'?); a similar range and 
mode as observed in this study. 
Estimated CPUE for all species 
combined in the longline fishery was 
0.0095 fish per hook-hour. Highest 
density CPUE (numbers of fish kept 
per 1,000 hook-hours) occurred in the 
eastern Gulf for red grouper and blueline 
tilefish, a similar distribution as reported 
by Ingram et al.'? In deeper waters of 
the western Gulf, yellowedge grouper, 
tilefish, and scamp had high CPUE den- 
sity values. For vertical line, the catch 
rate for all species was higher (0.0311 
fish per hook-hour) than observed in 
longline. Highest CPUE for red snapper 
occurred in the western Gulf, consistent 
with SEDAR.> Density CPUE values 


Ingram, W., M. Grace, L. Lombardi-Carlson, 
and T. Henwood. 2006. Catch rates, distribution 
and size/age composition of red grouper, Epi- 
nephelus morio, collected during NOAA Fish- 
eries Bottom Longline Surveys from the U.S. 
Gulf of Mexico. SEDAR-12-DW-05. Southeast 
Data Assessment and Review, South Atl. Fish. 
Manage. Counc., Charleston, SC (available at 
www.sefsc.noaa.gov/sedar/). 
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were higher and more dispersed in 
vertical line for other dominant species 
(vermilion snapper, red grouper, red 
porgy, and gag). 

As prescribed by NMFS’ National 
Bycatch Strategy addressing fishery 
bycatch on a national level, precision 
goals for bycatch estimates are defined 
in terms of CV estimates (NMFS, 2004). 
The precision of single species bycatch 
estimates is needed for population as- 
sessments; however, the reef fish fishery 
has bycatch from several stocks. In our 
study, CV estimates were low (0.1) for 
undersize target species, notably red 
grouper and red snapper. CV estimates 
for other species of commercial, rec- 
reational, and ecological importance, 
including several species of grouper 
and snapper, were relatively high and 
in some instances equal to 1.0. 

In terms of areas of high bycatch, 
management measures to reduce 
bycatch should consider targets that 
include changes in fishing behaviors 
relative to avoidance of high bycatch 
areas, modifications of gear to reduce 
bycatch, and cooperative efforts to 
close areas with high bycatch. As il- 
lustrated by Hot/Cold Spot Analysis!>, 
areas of highly significant rates of 
discards were identified. In longline, 
discard CPUE density was significantly 
higher in statistical areas 3 through 6. 
For vertical line, discard catch rates 
were significantly higher and concen- 
trated off Florida in statistical areas 5 
through 7. 

Prior to a mandatory observer pro- 
gram, self-reporting through logbook 
and discard supplementary data submis- 
sion were used to estimate sea turtle take 
projections in the reef fish fishery and 
formed the basis of biological opinions 
pursuant to formal consultation under 
Section 7 of the ESA (NMFS”°). Ob- 
servers documented twenty sea turtle 
interactions, notably in the bottom 


2°NMFS. 2005. Endangered Species Act-Section 
7 consultation on the continued authorization of 
reef fish fishing under the Gulf of Mexico Reef 
Fish Fishery Management Plan and Proposed 
Amendment 23. Biol. Opinion, 15 Feb., 115 
p. Southeast Reg. Off., Natl. Mar. Fish. Serv., 
NOAA, St. Petersburg. Fla. (available at http: 


longline component, during the study 
period (SEFSC'®), resulting in important 
implications for management. In Octo- 
ber 2009, a new biological opinion on 
the Gulf of Mexico reef fish fishery was 
completed with regulatory measures de- 
signed to minimize the impacts of future 
takes and monitor levels of incidental 
take (Fed. Regist.?'). 

Observer programs remain the most 
reliable means for monitoring fishery 
characteristics by not only providing 
insight on protected species interac- 
tions, but also for assessing quota and 
size restrictions, IFQ programs, CPUE, 
discard levels, gear effectiveness, and a 
wide array of other variables of interest 
to fishery managers, the fishing industry, 
academia, and the public. 
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Simulation of Tail Weight Distributions in Biological Year 1986-2006 
Landings of Brown Shrimp, Farfantepenaeus aztecus, 
from the Northern Gulf of Mexico Fishery 


CHARLES W. CAILLOUET Jr., RICK A. HART, and JAMES M. NANCE 


Introduction 


Size distribution within reported land- 
ings is an important aspect of northern 
Gulf of Mexico penaeid shrimp stock 
assessments. It reflects population char- 
acteristics such as numerical abundance 
of various sizes, age structure, and vital 
rates (e.g. recruitment, growth, and mor- 
tality), as well as effects of fishing, fish- 
ing power, fishing practices, sampling, 
size-grading, etc. (Kutkuhn, 1962; Neal, 
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1967; Rothschild and Brunenmeister, 
1984; Nance et al., 1994; Diop et al., 
2007; Caillouet et al., 2008; Nance et 
al., 2010; Parrack'; Nichols*). Age of 
shrimp cannot be determined directly 
(Parrack, 1979; Rothschild and Brunen- 
meister, 1984; Neal and Maris, 1985). 
Therefore, age structure of shrimp in 
reported landings has been determined 


'Parrack, M. L. 1981. Some aspects of brown 
shrimp exploitation in the northern Gulf of 
Mexico. Presented at the Workshop on the Sci- 
entific Basis for the Management of Penaeid 
Shrimp, Key West, Fla., Southeast Fisheries Sci- 
ence Center, National Marine Fisheries Service, 
NOAA, Miami, Fla. Unpubl. rep., 50 p. 
*Nichols, S. 1984. Updated assessments of 
brown, white and pink shrimp in the U.S. Gulf 
of Mexico. Presented at the Workshop on Stock 
Assessment, Miami, Fla., Southeast Fisheries 
Science Center, National Marine Fisheries Ser- 
vice, NOAA, Miami, Fla. Upubl. rep. 54 p. 


indirectly by estimating numbers of 
shrimp from pounds allocated to mar- 
keting size categories, and transforming 
size into age using growth curves (Neal, 
1967; Rothschild and Brunenmeister, 
1984; Nance et al., 1994; Parrack'!; 
Nichols?). 

Most but not all reported landings 
from northern Gulf of Mexico shrimp 
fisheries are size-graded. The usual mea- 
sure of shrimp size in archived landings 
data is count (C), the number of shrimp 
tails (abdomen or edible portion) per 
pound (0.4536 kg). Shrimp are mar- 
keted and landings reported in pounds 
within tail count categories. Statistically, 
these count categories are count class 
intervals or bins with upper and lower 
limits expressed in C. The upper and 
lower limits of most count class intervals 


ABSTRACT — Size distribution within re- 


ported landings is an important aspect of 


northern Gulf of Mexico penaeid shrimp 
stock assessments. It reflects shrimp popu- 
lation characteristics such as numerical 
abundance of various sizes, age structure, 
and vital rates (e.g. recruitment, growth, 
and mortality), as well as effects of fishing, 
fishing power, fishing practices, sampling, 
size-grading, etc. 

The usual measure of shrimp size in 
archived landings data is count (C) the 
number of shrimp tails (abdomen or edible 
portion) per pound (0.4536 kg). Shrimp 
are marketed and landings reported in 
pounds within tail count categories. Sta- 
tistically, these count categories are count 
class intervals or bins with upper and 
lower limits expressed in C. Count catego- 
ries vary in width, overlap, and frequency 
of occurrence within the landings. The 
upper and lower limits of most count class 
intervals can be transformed to lower and 
upper limits (respectively) of class intervals 


expressed in pounds per shrimp tail, w, the 
reciprocal of C (i.e. w = 1/C). 

Age based stock assessments have relied 
on various algorithms to estimate numbers 
of shrimp from pounds landed within count 
categories. These algorithms required un- 
derlying explicit or implicit assumptions 
about the distribution of C or w. However, 
no attempts were made to assess the actual 
distribution of C or w. Therefore, validity 
of the algorithms and assumptions could 
not be determined. When different algo- 
rithms were applied to landings within the 
same size categories, they produced differ- 
ent estimates of numbers of shrimp. 


This paper demonstrates a method of 


simulating the distribution of w in reported 
biological year landings of shrimp. We 
used, as examples, landings of brown 
shrimp, Farfantepenaeus aztecus, from the 
northern Gulf of Mexico fishery in biologi- 
cal years 1986-2006. Brown shrimp bio- 
logical year, T,, is defined as beginning on 
| May of the same calendar year as T; and 


ending on 30 April of the next calendar 
year, where subscript i is the place marker 


for biological year. Biological year land- 


ings encompass most if not all of the brown 
shrimp life cycle and life span. Simulated 
distributions of w reflect all factors influ- 
encing sizes of brown shrimp in the land- 
ings within a given biological year. Our 
method does not require a priort assump- 
tions about the parent distributions of 
w or C, and it takes into account the vari- 
ability in width, overlap, and frequency of 
occurrence of count categories within the 
landings. Simulated biological year distri- 
butions of w can be transformed to equiva- 
lent distributions of C. 

Our method may be useful in future test- 
ing of previously applied algorithms and 
development of new estimators based on 
statistical estimation theory and the under- 
lying distribution of w or C. We also exam- 
ine some applications of biological year 
distributions of w, and additional vari- 
ables derived from them. 
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can be transformed to lower and upper 
limits (respectively) of class intervals 
expressed in pounds per shrimp tail, w, 
the reciprocal of C (i.e. w = 1/C) 

Age based stock assessments have 
relied on various algorithms to estimate 
numbers of shrimp from pounds landed 
within count categories (e.g. Neal, 
1967; Rothschild and Brunenmeister, 
1984; Nance et al., 1994; Diop et al.. 
2007; Parrack'; Nichols”). These algo- 
rithms required underlying explicit or 
implicit assumptions about the distribu- 
tion of C or w. However, no attempts 
were made to assess the actual distribu- 
tions of C and w. Therefore, validity of 
the algorithms and assumptions could 
not be determined. When different algo- 
rithms were applied to landings within 
the same size categories (e.g. Parrack! 
vs. Nichols), they produced different 
estimates of numbers of shrimp (Cail- 
louet, 2003). 

Estimating numbers of shrimp from 
pounds landed within size categories is 
Statistically challenging for additional 
reasons. Some count categories rep- 
resenting the largest shrimp have an 
implied lower limit of zero (e.g. < I5 
count), and some representing the small- 
est shrimp have an implied upper limit 
of co (e.g. > 67 count). Neither zero nor 
eo can be transformed to real values of 
w. Count categories also exhibit con- 
siderable variability in width, overlap, 
and frequency of occurrence within 
the landings. Certain count categories 
dominate the landings, reflecting what 
are referred to as standard count catego- 
ries: <15, 15-20, 21-25, 26-30, 31-40, 
41-50, 51—67., and > 67 count (Caillouet 
etal., 2008). 

This paper demonstrates a method 
of simulating the distribution of w in 
reported biological year landings of 
shrimp, as a basis for further investiga- 
tion and evaluation of previously used 
algorithms and development of new 
ones. We used, as examples, landings 
of brown shrimp, Farfantepenaeus 
aztecus, trom the northern Gulf of 
Mexico fishery in biological years 
1986-2006. Neal (1967) defined brown 
shrimp biological year, T,,as beginning 
1 May of the same calendar year as T. 


I 


and ending 30 April of the next calen- 


Table 1.—Symbols and descriptions of variables used in analyses of biological year reported landings of brown 
shrimp from the northern Gulf of Mexico fishery. These apply only to size-graded landings in legitimate count 
categories; i.e. data selected by filtering, editing, and removing residual outliers from archived landings data. 


Symbols 


Descriptions of variables 


A biological year, from 1 May of a given calendar year through 30 April of the next calendar year, where /= 
0 20 is the place marker for biological years 1986-2006 


Cc. the /” lower limit of a legitimate count (number per pound) category in landings data from the ?” biological 


year, where / = 0 m, 


the total number of C, in landings data from the 7” biological year 
the /” upper limit of a pounds per shrimp tail category, where w, = /C;, in landings data from the /” 


biological year 


the 7” cumulative proportion of pounds landed at w, in ?” biological year 


the /” weighting factor for the P, and w, data pairs in the ?” biological year. This weighting factor, q,, is 
the sum of observations over all count categories having C, as their lower limit (or w, as their upper limit), 
regardless of the recorded upper limits of these count categories 


the k” simulated value of weight per shrimp tail, where 0.005155 Ib < w’, < 0.111111 Ib, k=0,..., 999, and 


the interval between the w’, is 0.000106 


the «” cumulative proportion of pounds landed at w’, in the ”” biological year, which is simulated from the 
modified Richards function fitted to P, on w, in the ? biological year 


the parameter, estimated from the modified Richard's function fitted to P, on w, in the /” biological year, 
which allows the w’, at which P’, = Pmax,/2 to vary among biological years 


the parameter, estimated from the modified Richard's function fitted to P, on w, in the 7” biological year, 
which represents the maximum intrinsic rate of increase in P’, per unit w‘, at the inflection point of the curve 
the parameter, estimated from the modified Richard's function fitted to Pon w, in the # biological year 
which allows the sigmoid shape of the curve to vary (symmetrical or asymmetrical) among biological years 
the k” simulated proportion of pounds landed at w’, in the 7” biological year 


the ?” biological year yield, which includes pounds of brown shrimp tails landed in legitimate count 
categories and in the unknown size category combined 


the k” simulated number of shrimp tails at w’,, where 0.005155 Ib < w’, < 0.111111 Ib, in the #* biological 


year 


the simulated total number of shrimp tails landed in the ?” biological year 


the simulated pounds per shrimp tail at which half of Y, is harvested in the ?” biological year 
the simulated mean count of brown shrimp in the landings from the /” biological year 


the simulated mean pounds per shrimp tail of brown shrimp in the landings from the i" biological year 


dar year, where subscript / is the place 
marker for biological year (Table 1). 
Most landings in the i” biological year 
are assumed to be produced from co- 
horts recruited to the fishery within that 
same biological year. In other words, 
a biological year encompasses most of 
the cycle and life span of brown shrimp 
within this intensive fishery. 

Our approach does not require a 
priori assumptions about the parent 
distributions of w or C, and it takes into 
account the variability in width, overlap, 
and frequency of occurrence of count 
categories within the landings. Simu- 
lated biological year distributions of w 
can easily be transformed to equivalent 
distributions of C. Our method may be 
useful in future testing of previously 
applied algorithms and development 
of new estimators based on statistical 
estimation theory and the underlying 
distribution of w or C. We also examine 
some applications of biological year dis- 
tributions of w and additional variables 
derived from them. 


Materials and Methods 
Fishery 


The brown shrimp fishery of the 
northern Gulf of Mexico is bounded 
by statistical subareas 10-2), and com- 
prises inshore (estuarine) and offshore 
(Gulf of Mexico) territorial waters of 
Texas, Louisiana, Mississippi, Ala- 
bama, and a portion of Northwestern 
Florida, as well as adjoining Federal 
waters landward of the 50 fm depth 
contour within the U.S. Exclusive 
Economic Zone (EEZ) (Fig. 1). Brown 
shrimp produce annual crops (Neal and 
Maris, 1985), with recruitment to the 
fishery occurring in May—July (Roth- 
schild and Brunenmeister, 1984). 
though life span is 20-27 mo (Baxter, 
1971), most brown shrimp are har- 
vested within 6 mo of age.’ Neal (1967) 


‘Fishery Management Plan for the Shrimp Fish- 
ery of the Gulf of Mexico, United States Waters. 
Gulf of Mexico Fishery Management Council, 
Tampa, Fla., Nov. 1981 (online at http://www. 
gulfcouncil.org). 
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Figure 1.— Shrimp Statistical Subareas 10-21, encompassing the brown shrimp fishery within inshore (estuarine) and offshore (Gulf 
of Mexico) state territorial waters, and part (within the 50 fm depth contour) of the adjoining Federal EEZ in the northern Gulf of 


Mexico. 


conducted virtual population analyses 
of brown shrimp in statistical subareas 
18 and 19 (Fig. 1), and found that esti- 
mated numbers of brown shrimp in re- 
ported landings during biological year 
1964 represented 97.7% of the total 
virtual population over a 17-mo period. 
This finding indicated that only 2.3% 
(by number) of the shrimp recruited as 
new cohorts in biological year 1964 
contributed to the landings in biological 
year 1965. If shrimp landed in a given 
biological year within our time series 
(1986-2006) included survivors from 
cohorts recruited in preceding biologi- 
cal years, this could have affected our 
biological year simulations of w and 
other variables derived from them. 
However, such a carryover would be 
small, because it would involve only 
the larger sizes of shrimp which are 


73(2) 


lowest in pounds and fewest in numbers 
within the landings. 


Landings Data 


Brown shrimp landings data are 
archived by the National Marine 
Fisheries Service (NMFS) Galveston 
Laboratory, Texas. Statistically, re- 
ported landings are fishery-dependent 
samples taken without replacement 
from the brown shrimp population. 
They are multitudinous but have limi- 
tations (Kutkuhn, 1962: Snow, 1969; 
Prytherch, 1980; Parrack'; Nichols’; 


Poffenberger*) which may bias not 


*Poffenberger, J. R. 1991. An overview of the 
data collection procedures for the shrimp fish- 
eries in the Gulf of Mexico, National Marine 
Fisheries Service, Southeast Fisheries Science 
Center, Miami, Fla. (online at http://www.setsc. 


noaa.gov/gssprogram.jsp). 


only our simulated distributions of w 
and additional variables derived from 
them, but also may have biased previ- 
ous estimates of numbers of shrimp 
from pounds landed within count cat- 
egories. Not all brown shrimp that are 
caught are landed, and not all that are 
landed are reported (Kutkuhn, 1962; 
Berry and Benton, 1969; Baxter, 1973; 
Snow, 1969; Prytherch, 1980, Nance 
et al., 1991; Caillouet et al., 2008; 
Poffenberger*). Nonreported catch 
includes shrimp marketed directly to 
consumers, marketed as fishing bait 
(not all, but some), discarded for vari- 
ous reasons, Kept for personal use by 
shrimpers, or otherwise not reported. 
Thus, reported landings are less than 
the actual catches, and also represent 
incomplete samples of the actual land- 
ings (Caillouet et al., 2008). 
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Reported shrimp landings data are 
recorded by calendar year, month, 
statistical subarea (Fig. 1), depth zone, 
shrimping trip, and count category or 
unknown size category, along with other 
information (Kutkuhn, 1962; Snow, 
1969; Prytherch, 1980; Poffenberger*). 
We treated the unknown size category 
as a catch-all category. In selecting re- 
cords for a working file of size-graded 
landings data for our simulations, we 
excluded all landings originally reported 
in the unknown size category, as well 
as landings added to the unknown size 
category after we judged their count cat- 
egories to be outliers (see Data Selection 
and Preparation below). The resultant 
unknown size category contained land- 
ings that were: 


1) not size-graded, 


2) size-graded incorrectly or size 
limits not recorded, 

3) not assigned to a count category 
for other reasons (e.g. pieces of 
shrimp tails), or 

4) size-graded but reported in count 


categories we judged to be outliers. 


Previous investigators (e.g. Roths- 
child and Brunenmeister, 1984; Par- 
rack'; Nichols*) also excluded certain 
landings from their analyses for vari- 
ous reasons. Two methods of grading 
shrimp, box-grading and machine 
grading, were described by Kutkuhn 
(1962), Snow (1969), Prytherch (1980), 
and Poffenberger*. Differences between 
these grading methods and variations 
in their relative contributions to size- 
graded landings over time may have 
biased our simulated distributions of w 
and variables derived from them, but 
they may also have biased previous 
estimates of numbers of shrimp within 
count categories. 


Data Selection and Preparation 


Our final working file contained ar- 
chived landings records selected from 
biological years 1986-2006, but only 
those we considered to have legitimate 
count class limits. We initially consulted 
NMFS port agents (who collect landings 
data) to obtain their opinions about the 
true range in size of brown shrimp tails 


in the landings. It was agreed that the 
maximum C (smallest shrimp) for brown 
shrimp in the landings was around 250 
tails per pound (equivalent to w = 0.004 
Ib, or 1.8 g), and minimum C (larg- 
est shrimp) around 9 tails per pound 
(equivalent to w =0.111 Ib, or =50.3 g). 

Preparation of the working file in- 
volved filtering and editing a copy of 
archived data from biological years 
1986-2006 as follows: 


1) If a record was originally coded 
as belonging to the unknown cat- 
egory, it was excluded. 

If an upper or lower limit of a 
count category was not recorded 
(i.e. left blank), the record was 
excluded. 

If a recorded lower limit exceeded 
the recorded upper limit of a count 
category, the limits were assumed 
to have been inadvertently trans- 
posed at data entry, and the record 
was retained in the working file 
after being recoded by interchang- 
ing its count category limits. 

If recorded upper and lower limits 
of a count category were both C = 
0, the record was excluded. 

If the recorded upper limit of a 
count category was 0 < C <9, both 
the lower and upper limits were 
recoded as C = 9, and the record 
was retained in the working file. 
If only the recorded lower limit of 
acount category fell within C <9, 
but the recorded upper limit was > 
9, the lower limit was recoded as 
C =9, and the record was retained 
in the working file. 

If the recorded lower and upper 
limits of a count category were C 
> 250, the record was excluded. 
If the recorded upper limit of a 
count category was C > 250, but 
the recorded lower limit was C < 
250, the upper limit was recoded 
as C = 250, and the record was 
retained in the working file. 

All other archived records were 
retained in the working file. 
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We then performed statistical analy- 
ses of the working file to identify and 
remove records having count class 
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Figure 2.— Upper vs. lower limits of 
brown shrimp count categories in fil- 
tered and edited landings in biologi- 
cal year 2006; (A) before residual 
outlier records were removed and 
(B) after residual outlier records 
were removed. Lines were fitted by 
weighted linear regression, where 
the weighting factor was number of 
shrimping trips associated with each 
unique count category. 


limits we judged to be outliers. For each 
biological year, we used SYSTAT? to fit 
preliminary weighted linear regressions 
of upper limits on lower limits of the 
count categories, where the weighting 
factor was the number of observations 
(i.e. shrimping trips) associated with 
each unique count category (i.e. unique 
combination of upper and lower limits). 
Figure 2A is an example of a preliminary 
regression and data plot for biological 


5Mention of trade names or commercial firms 
does not imply endorsement by the National 
Marine Fisheries Service, NOAA. 
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year 2006. Statistical weighting by 
number of shrimping trips was our way 
of dealing with variability in frequency 
of occurrence of count categories in the 
working file. Records removed from the 
working file by filtering, editing, and 
identification of residual outlier count 
categories represented a higher percent- 
age of observations than percentage 
of pounds landed (Table 2): i.e. they 
contained relatively low pounds per 
observation. 

We fitted final weighted linear regres- 
sions of upper limits on lower limits 
within the final working file for each 
biological year (Table 3). The weight- 
ing factor for these regressions was the 
number of observations (i.e. shrimping 
trips) associated with each unique count 
category remaining in the final working 
file. These final regressions character- 
ized the relationship between legitimate 
count category upper and lower limits 
for each biological year. Figure 2B is an 
example final regression and data plot 
for biological year 2006. Slopes and 
intercepts of the final linear regressions 
(Table 3) for each biological year were 
examined for trends, using polynomial 
regression. Coded biological year (T; 
— 1996) was substituted for 7, in these 
polynomial regressions, to avoid prob- 
lems that otherwise might have been 
caused by correlations among powers 
of T, (Sokal and Rohlf, 2000). 


Aggregation and Cumulation 
of Landings 


Landings from the final working file 
were aggregated (summed) by biologi- 
cal year and count category lower limits, 
C;,, where j is the place marker for the 
C;, within a biological year; j = 0...., 
m;, Where m, is the total number of C;; 
in each biological year (Table 1). Upper 
limits of count categories (equivalent to 
lower limits of class intervals of w) were 
ignored. Because the number of unique 
count categories in the final working 
file varied among biological years, 
the number of C;; also varied among 
biological years, as did m;. Summing 
the landings by biological year and C;; 
produced a subset of data with much 
lower spatial-temporal resolution than 
that of more detailed data sets used 


Table 2.—Number of observations (shrimping trips) 
and pounds landed in the NMFS-archived records, 
compared to those remaining after filtering, editing, 
and removal of residual outlier count categories, for 
brown shrimp landings in the northern Gulf of Mexico 
fishery in biological years 1986-2006. 


Observations 
Records (shrimping trips) Pounds (tails) 
Archived 2,425,373 1,682,806,769 
100.0% 100.0% 
After filtering and 2,319,554 1,668,305,100 
editing 95.64% 99.14% 
Aiter residual outlier 2,308,674 1,664 ,449,467 
removal 95.19% 98.91% 


in previous, bottom-up approaches to 
estimating numbers of shrimp within 
count categories (Neal, 1967; Roths- 
child and Brunenmeister, 1984; Nance 
et al., 1994; Diop et al., 2007; Parrack'; 
Nichols’). 

Biological year summations of land- 
ings combined all spatial-temporal influ- 
ences (statistical subarea, depth zone, 
and month) on size of brown shrimp in 
the landings. These influences included 
sex ratio, recruitment, growth, mortality, 
fishing effort, fishing power of shrimp 
trawlers, experience of captains and 
crews, gear selectivity, discarding, data 
collection procedures, grading methods, 
and possibly other factors that affect 
count category landings within a bio- 
logical year. Spatial influences were col- 
lapsed to the level of the entire fishery, 
and temporal influences to the level of 
biological years. Summation of landings 
by C,, combined landings within count 
categories having C;, as their lower limit. 
The simple hypothetical example below 
depicts this process: 


Count Pounds 
Category Observations landed 
9-12 2 500 
9-15 1 200 
9-20 l 40 
Total 6 1,740 


The sum of observations over all 
count categories having C,, as their 
lower limit became the weighting factor, 
q;;, for each C;, and the sum of pounds 
associated with it. In the hypothetical 
example above, C;,=9,q,, = 6, and both 
are associated with 1,740 lb landed. 


Table 3.—Final weighted linear regressions of upper 
(U) on lower (L) limits of count categories in brown 
shrimp landings data selected by filtering, editing, and 
removal of residual outliers from the NMFS-archived 
landings data. The weighting factor was the number 
of shrimping trips associated with each unique count 
category (i.¢. unique U and L data pair) in the landings 
data selected from each biological year. Sample size 
was the sum of these weighting factors for each bio- 
logical year (see Fig. 2). 


Biological 


Sample Adjusted 
year, T, Intercept, Slope, size r2 
1986 —0.0609420 1.172878 141,523 0.988 
1987 -0.7467180 1.189633 159,010 0.988 
1988 0.4795405 1.160103 158,733 0.992 
1989 0.0479237 1.171595 147,315 0.992 
1990 0.6609263 1157006 137,647 0.993 
1991 -0.1564869 1.180132 122065 0.992 
1992 0.1879965 1.169404 117,633 0.991 
1993 -0.5079871 1.187639 105,907 0.989 
1994 0.0533226 1.173414 111,968 0.992 
1995 -0.1264397 1.179191 102643 0.993 
1996 —0.7873293 1.195364 97,111 0.989 
1997 —1.4422680 1.210573 98,415 0.987 
1998 -1.0609130 1.199557 91,378 0.988 
1999 —1.2453040 1.204808 92,638 0.985 
2000 -0.3466840 1.179789 95,775 0.990 
2001 0.1584804 1.169094 89,022 0.992 
2002 —0.2696291 1.187891 122,160 0.992 
2003 -1.0236740 1.206445 103,013 0.993 
2004 -1.1202140 1.208740 82,006 0.993 
2005 —0.3479283 1.192750 69,662 0.994 
2006 -0.9036610 1.208740 63,050 0.995 


Examples of variation in C;, and q;, for 
biological years 1986, 1996, and 2006 
are shown in Figure 3. Dominant C; 
were conspicuous as indicated by their 
q;;, and many were identical or close to 
the C;, of standard count categories, as 
expected. 

Within each biological year, the 
pounds associated with C;, were cumu- 
lated over the observed range of C,, 
from the highest to the lowest C;, (i.e. 
from the smallest to largest shrimp tails). 
These cumulative pounds were then 
converted to proporticns of cumulative 
pounds landed, P;, (Table 1), from the 
highest to the lowest C;,. Figure 4A is an 
example of the stair-stepped relationship 
between P,, and C,, for biological year 
2006, and Figure 4B is the equivalent 
stair-stepped relationship between P,, 
and where 1/C;,. 

Modified Richards Function 

We searched for an asymptotic, asym- 
metrical sigmoid regression model to 
convert the stair-stepped relationship 
between P;, and w,, to a smooth curve 
for each biological year. The regression 
model we chose was a simplified form of 
the Richards function (Richards, 1959): 
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Figure 3.— Weighting factors, q,, (i.e. shrimping trips) vs. legiti- 
mate brown shrimp count category lower limits, C,,. for biologi- 
cal years (A) 1986, (B) 1996, and (C) 2006. Dominant C j are 
marked by the numbers above vertical bars representing their qj;- 


P = Pmax(l-e*”") 


where, P 


(1) 


is the cumulative proportion 
of pounds landed at w, 


» is shrimp tail weight in 


pounds, over the observed 
range from minimum to max- 
imum w, 

is the upper asymptote, 

is the parameter which allows 
wat which P= Pmax/2 to vary, 
is the parameter which repre- 
sents the maximum intrinsic 
rate of increase in P per unit 
w, which occurs at the inflec- 
tion point on the curve, 


* is the parameter that allows the sig- 
moid shape of the curve to vary (sym- 
metrical or asymmetrical), and 
is the base of natural logarithms. 


Because we constrained Pmax to 
equal | in fitting all the regressions, Eq. 
(1) was simplified into the following 
regression model: 


P=(l-e*") 


For each biological year, we used 
GraphPad Prism (version 5.02) to fit 
Eq. (2) to P;,on w,, by weighted non- 
linear regression, where the weighting 


Figure 4.— Biological year 2006 (A) cumu- 
lative proportion of pounds landed, P59), 
vs. lower limit. C,,;, of count categories of 
filtered and edited brown shrimp landings 
from which residual outlier records were 
removed and (B) relationship between P 


20) 
and where 1/C 


20;" 


factor was q;,,. In this way, parameters 
a;, b,, and c,; (Table 1) were estimated 
for each biological year (Table 4). The 
lower case parameter c, should not be 
confused with the upper case count C;,. 
We tried fitting a number of other asym- 
metrical sigmoid functions available in 
GraphPad Prism, but Eq. (2) was the 
best fitting of those we examined. While 
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we recognize that additional curve fit- 
ting methods and models could have 
been tested, Eq. (2) was adequate for 
purposes of demonstrating our simula- 
tion approach. By fitting Eq. (2), we 
smoothed the relationship between P,, 
on w;,, and obtained an equation repre- 
senting this relationship for each biolog- 
ical year (Table 4). We also calculated 
the adjusted r? as an approximation 
of how well Eq. (2) fit the data points 
for each biological year (Table 4), but 
recognize it is not strictly applicable to 
nonlinear regression. 


Simulating Biological Year 
Distribution of Tail Weight 


The next step toward simulating the 
distribution of w was to generate a new 
set of data pairs for each biological 
year, using the fitted equations in Table 
4. First, we generated equally spaced 
values of w’, (Table 1), from a minimum, 
w’, (= 0.005155 Ib), to a maximum, 
W’ ggg (=0.111111 Ib), where the k"" place 
marker for the w’, was k = 0, ... , 999 
(Table 1). The increment, g, between the 
w’, was then calculated as 


( Wooo )/999 
= 0.000106 Ib. 


The w’, were generated by 
wi =k(g)+wy. 


We then generated values of P’,, for 
each w’, for each biological year, using 
the following equation and estimates of 
parameters a; , b;, and c; from Table 4: 


(3) 


Three reasons for applying w% = 
0.005155 Ib (derived from 1/194) as 
the minimum shrimp tail weight for all 
biological year simulations were: 


1) The lowest maximum C;, observed 
(in the working file) among all 
biological years was 194 count, 
the reciprocal of the highest mini- 
mum 

Imaginary numbers were gener- 
ated by Eq. (3) for the minimum 
P’,,in some biological years when 
the actual minimum w;; observed 


2) 


Table 4.—Biological year yield (Y), parameter estimates, and other statistics for weighted nonlinear regressions 
(modified Richards function, Eq. (2)) of cumulative proportions of pounds landed, P,, on pounds per shrimp tail, 
w;, in brown shrimp landings data selected by filtering, editing, and removal of residual outliers from the NMFS- 
archived landings data. The weighting factor was the number of shrimping trips, Qip associated with each data 


pair, P,, and w,, in the selected landings data. For a given biological year, the number of data points analyzed (total 


sample size) was the sum of these weighting factors, m 


0 
Total sample 
size 
Estimated parameters m, 

Biological Yield Y,, Ya, Adjusted 
year, T, pounds a, b, j=0 
1986 94,738,424 0.2770842 53.75477 1.016937 141,523 0.994 
1987 89,394,421 0.3177634 61.64658 1.074417 159,010 0.997 
1988 79,859,436 0.2713897 57.15293 1.286208 158,733 0.997 
1989 94,170,525 0.2802385 58.89570 1.243767 147,315 0.996 
1990 105,121,282 0.2865627 55.59358 0.968642 137,647 0.992 
1991 85,602,708 0.1627544 47.08183 1.095961 122,065 0.993 
1992 68,425,417 0.2294646 55.76027 1.150789 117,633 0.995 
1993 66,431,237 0.2427682 55.10823 0.865503 105,907 0.989 
1994 67,049,354 0.2126820 51.46945 1.107677 111,968 0.996 
1995 75,859,021 0.2123855 48.21137 0.829590 102,643 0.991 
1996 73,500,416 0.2459783 55.83692 0.888528 97,111 0.991 
1997 65,389,618 0.2837078 55.32308 0.761172 98,415 0.994 
1998 80,514,861 0.2723718 61.82822 0.975136 91,378 0.992 
1999 81,035,496 0.2308989 56.10879 0.836558 92,638 0.987 
2000 94,463,851 0.3038908 59.25881 1.084649 95,775 0.995 
2001 87,660,251 0.3287329 74.62214 1.352838 89,022 0.987 
2002 73,180,653 0.3917993 81.88587 1.447248 122,160 0.988 
2003 82,309,001 0.3194503 79.86258 1.376817 103,013 0.986 
2004 74,233,767 0.2973424 57.98183 0.884165 82,006 0.981 
2005 58,819,403 0.2349768 56.86499 1.169126 69,662 0.981 
2006 85,047,627 0.0818478 51.68214 1.640370 63,050 0.991 


in those years was applied (this 
probably was due in part to the 
fact that Eq. (3) did not fit the data 
points representing very small 
shrimp tails closely in those years). 
It was consistent to constrain w’, 
to be the same for all biological 
years. 


Ww 
— 


The first derivative of Eq. (3), 6P’,/8’,,,, 
was 


For each biological year, we used Eq. 
(4) to generate first derivatives for each 
w’,. To transform these first derivatives 
(Eq. (4)) into p’, (Table 1), which was 
the proportion of pounds landed at w’, 
for each biological year, we divided 
them by the sum of all first derivatives 
over the range in w’,, for each biological 
year. This sum was calculated as 


999 


> 


k=0 


In other words, for each biological year, 
p’,, at each w’, was calculated as 


= (SP, /dw7). 


k=0 


Biological year yield, Y,, encompassed 
all landings within a biological year, in- 
cluding those retained in our final work- 
ing file as well as those that had been 
excluded from it. For each biological 
year, number of shrimp tails, /’,, (Table 
1), at each w’, was calculated by 


- (5) 


Equation 5, describing the relationship 
between /’, and w’,, is the simulated dis- 
tribution of w for the i” biological year. 

We would have been able to exclude 
some steps in our simulation sequence 
had the final working file represented 
total reported landings from each 
biological year (i.e. Y,). However, 
the final working file was a subset of 
size-graded landings selected from the 
archived landings, and it did not con- 
tain landings we excluded (i.e. those 
relegated to the unknown category), 
whereas Y, contained all landings for 
each biological year. Therefore, Eq. (5) 
applied the subset of proportions p’,, to 


= 
a 
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the total yield Y, to estimate /’, for each 
biological year. 

We recognize that relative distribu- 
tions of w for each biological year, and 
their corresponding cumulative relative 
distributions, also could have been de- 
rived from our simulated distributions 
of w. They might be of interest in some 
applications of our approach, but they 
were not essential to the purpose of our 
paper. They can easily be calculated 
from the information provided in this 
paper. However, the concept of cumu- 
lative relative distribution of w in bio- 
logical year landings of brown shrimp 
is important in that it would estimate 
the probability of occurrence of tail 
weight <w; i.e. it would be an approxi- 
mation of the cumulative distribution 
function (CDF) for w. This is the major 
part of the explanation of why we chose 
lower limits, C;, (equivalent to upper 
limits of w..), for aggregating and cumu- 
lating landings, and then transformed 
C;; to w;; in preparation for fitting Eq. 
(2 ). Because a simulated distribution 
of w can be used to calculate the rela- 
tive distribution of w and cumulative 
relative distribution of w, it is relevant 
to future testing of past algorithms and 
development of new ones to estimate 
numbers of shrimp from pounds landed 
within class intervals of w or C in the 
landings. Although we excluded certain 
landings (unknown size category) and 
ignored upper limits of legitimate count 
categories in simulating biological year 
p’,. our simulations of f’, included all 
biological year landings (Y,); i.e. all 
biological year landings contributed to 
simulation of biological year distribu- 
tions of w. 


Biological Year 

Total Number of Shrimp 
Tails (N;), Mean C,, 

and Mean w, 


The total number of shrimp tails, N;. 
in the landings from a biological year T, 
was simulated by 


(6) 
k=0 


Crude estimates of biological year 
mean count (N,/Y,) and its equivalent 


mean tail weight (Y/N;) were calcu- 
lated. We examined trends in both of 
these means via polynomial regression, 
where coded years (T, — 1996) were 
substituted for T.. 


Tail Weight at Half of Y, 


Given that a fitted equation — 
ing the relationship between P;, and w;, 
was available for each biological year 
(Table 4), we estimated tail weight, w50,, 
at which half of the annual yield, Y,/2, 
was harvested in each biological year 
(note that when Pmax is constrained 
to equal 1, w50; = Pmax/2 = 0.5). Each 
equation (Table 4) was solved for w50; 
as follows: 


w50, 


This statistic is similar in concept to 
LDSO, the estimated lethal dose (con- 
centration) of a toxic substance at 
which 50% mortality occurs in exposed 
subjects. In our application, it is a po- 
tentially useful index of the relationship 
between brown shrimp size and yield 
(see Caillouet et al., 2008). We exam- 
ined w50. via polynomial regression, 
where coded years (7, — 1996) were 
substituted for T.. 


Results 


Polynomial Regressions 


We recognize that polynomial re- 
gression is an empirical approach to 
fitting a curve to a time series of data, 
and that the resulting polynomial terms 
have no structural meaning (Sokal and 
Rohlf, 2000). We applied it only to 
detect possible trends in the variables 
we simulated, and to demonstrate 
possible applications of our simulated 
distributions of w. Obviously, many 
other curve fitting approaches could 
have been used to examine the time 
series for each variable. Causes and 
effects within this brown shrimp fish- 
ery could have influenced the detected 
polynomial trends, despite variability 
(deviations from regression) caused 
by fluctuations in annual recruitment 
and other factors which are typical in 
shrimp populations (Caillouet et al., 
2008). 


Weighted Linear Regressions 
of Upper vs. Lower 
Limits of Count Categories 


Data plots and preliminary weighted 
linear regressions of upper on lower 
limits of unique count categories in 
each biological year (see the example 
for the year 2006 in Fig. 2A) showed 
that count category outliers remained 
in the data after filtering and editing. In 
the year 2006 example, the outliers were 
concentrated near the minimum lower 
limit count of 9 (largest shrimp), which 
elevated the intercept of the fitted line in 
Figure 2A as compared to the intercept 
of the fitted line in Figure 2B, in which 
residual outliers had been removed. The 
unusually wide class intervals of outlier 
count categories could lead to serious 
biases in estimating numbers of shrimp 
within such count categories. Also, we 
emphasize that each data pair (upper 
and lower limits) was weighted, so the 
actual numbers of residual outliers are 
much higher than the number of data 
points representing outliers in Figure 
2A (Table 2). 

As expected, final weighted linear 
regressions of upper limits on lower 
limits of count categories were close 
fitting in all biological years as shown 
by high adjusted r? (Table 3, Fig. 2B). 
These final regressions characterized 
the relationship between upper and 
lower limits of what we considered to 
be legitimate count categories in each 
biological year. All slopes of these final 
regressions were slightly greater than 
1 (Table 3), indicating that count class 
intervals in the working file widened as 
their lower limits increased. Trends in 
slopes and intercepts of these regres- 
sions are shown in Figures 5A and 5B, 
respectively. 


Biological Year m; 
and Weighted 
Regressions of 


The biological year total number, 

, of C;, exhibited a concave quadratic 
(parabolic) trend (Fig. 6); m; dropped 
from 77 in 1986 to 35 in 1995, then 
increased but not to its earlier highest 
level. This trend in m; reflected changes 
in the total number of legitimate 
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intercept, = -0.4075 - 0.04810(T,- 1996) 
adjusted r? = 0.219 


Intercept, 


Slope, = 1.186 + 0.001723(T, - 1996) 
adjusted r* = 0.387 
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12 3 0 4 8 12 
T,- 1996 


Figure 5.—Linear trends in (A) 
intercepts and (B) slopes of final 
weighted linear regressions of upper 
limits on lower limits of legitimate 
brown shrimp count categories, over 
coded biological years (7, — 1996). 


count categories over the biological 
years. However, the total number of 
count categories in biological year T, 
exceeds m,, because upper limits of 
count categories were ignored in our 
simulations; i.e. landings at the count 
category level were combined at the 
count category lower limit level, C;; 
(see Aggregation and Cumulation of 
Landings). Wide variation in biological 
year numbers of count categories and 
the consequential quadratic trend in 
m, (Fig. 6) are interesting and worthy 
of further investigation. They could 
reflect changes in size-related market- 
ing strategies, recruitment, and perhaps 
other influences on choices of count 
categories in the landings. 

Weighted nonlinear regressions of 
P;, on w, for all biological years were 
close fitting, as indicated by very high 
adjusted r7 (Table 4). Over all biological 
years, adjusted r? equaled or exceeded 
0.981. Examples of plotted data points 
P;, vs. w;, and fitted curves for 1986, 
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Im, = 42.26 - 0.07792(T, - 1996) + 0.2421(T, - 1996)" 
adjusted r° = 0.523 


70- 


0 
T,- 1996 


Figure 6.—Quadratic trend in total 
number, m;, of count category lower 
limits, C;, for filtered and edited 
brown shrimp landings from which 
residual outlier count categories were 
removed, over coded biological years 
(T, — 1996). 


1996, and 2006 are shown in Figures 
7A-C, respectively. Inflection points of 
the regressions were far to the lower left 
in such plots (Fig. 7A—C), suggesting 
that brown shrimp were fully recruited 
to the landings at very small sizes, which 
is a very important finding. 
The total sample size, 


m, 


j=0 


(Fig. 8), for each biological year re- 
gression (Table 4), and the adjusted r7 
(Fig. 9) for these regressions, declined 
over biological years. In other words, 
adjusted r7 and total sample size were 
dependent, as expected (Fig. 10); i.e. 
the larger the sample size the higher the 
adjusted r2. We emphasize that the total 
sample size (Fig. 8) used in fitting the 
regressions of P,,on w;, for each biologi- 
cal year was less than the actual number 
of shrimping trips in the archived data 
for each biological year, because land- 
ings from some trips were initially in the 
unknown category or later placed there 
by filtering, editing, and outlier removal 
from the working file. Therefore, the 
data points and trend in Figure 8 should 
not be taken to represent total shrimping 
trips in the biological years. 

As is common in fitting models 
containing more than one parameter, 
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Figure 7.— Weighted nonlinear re- 
gressions (modified Richards func- 
tion) fitted to cumulative proportion 
of brown shrimp landings, P;,, vs. 
pounds per shrimp tail, w,,, for bio- 
logical years (A) 1986, (B) 1996, and 
(C) 2006. 


the parameter estimates often are not 
independent (i.e. orthogonal). Graph- 
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Figure 8.— Linear trend in an- 
nual shrimping trips, 


m 


j=0 
over coded biological years 
(T, — 1996), for brown shrimp 
landings in legitimate count 
categories. 


Pad Prism provided estimates of de- 
pendency of estimated parameters a,, 
b,, and c; within each biological year 
regression (dependency = | represents 
complete dependency, and dependency 
= 0 indicates orthogonality). Over bio- 
logical years, dependency was 0.865-— 
0.985 for parameter a,, 0.968-0.981 
for parameter b,, and 0.978-0.994 for 
parameter c;. Not only did all these 
parameters show strong dependency 
within each biological year regres- 
sion, but they also appeared related to 
each other over biological years (Fig. 
11A-C). 


Simulated Distributions of w 

Example distributions of w’, for 
biological years 1986, 1996, and 2006 
are shown in Figures 12A-C. All were 
strongly skewed to the right. Their 
most striking feature was their like- 
ness to negative exponential curves. 
Therefore, we plotted them in the form 
of In(f’,,) vs. w’, for all biological years 
(Fig. 13). Straight lines for In(f’,,) vs. 
w’, would have indicated that these 
simulated distributions of w followed 
a negative exponential pattern, once 
full recruitment to the landings was 
reached at very small sizes (Fig. 13). 
Only slight concavity was evident in 
all the curves. 


adjusted r* "0.9911 - 0.0005675(T, - 1996 
adjusted ~ = 0.530 


Adjusted r;” 


0 
T, - 1996 


Figure 9.—Linear trend in adjusted 
r? for weighted nonlinear regressions 
(modified Richards function) fitted 
to cumulative proportion of brown 
shrimp landings, P,,. vs. pounds per 
shrimp tail, wij, over coded biological 
years (T, — 1996). 
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Figure 10.—Linear relationship be- 
tween adjusted 7 and annual shrimp- 
ing trips, 
qi. 
j=0 
for weighted nonlinear regressions 
(modified Richards function) fitted 
to cumulative proportion of brown 
shrimp landings, P;,, vs. pounds per 
shrimp tail, w,.. 


Biological Year Total Number 
of Shrimp Tails and Yield 


Interestingly, although the biological 
year total number of shrimp tails, N; 
(Fig. 14), and yield, Y, (Fig. 15), showed 
hints of declines, they exhibited no 
significant trends over biological years, 
because of wide year to year variation. A 
close linear relationship between N; and 


b, = 61.36 - 171.84, + 587.0 


adjusted = 0.731 


0.3948 + 0.01195b; 


adjusted Fa = 0.187 


= 2.396 - 11.694, + 24.314” 


adjusted Fad = 0.488 


Figure 11.—(A) Quadratic relation- 
ship between parameters b; and q; . 
(B) linear relationship between par- 
ameters c, and b., and (C) quadratic 
relationship between parameters c, 
and a,, for weighted nonlinear re- 
gressions (modified Richards func- 
tion) fitted to cumulative proportion 
of brown shrimp landings, P;,, vs. 
pounds per shrimp tail,w,. 
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Y, (Fig. 16) was expected; i.e. the more 
pounds landed the greater the number 
of shrimp tails in the landings, and vice 
versa. However, biological year mean 
count, N,/Y, (Fig. 17A) was not constant, 
because simulated distributions of w 
and Y, were not constant over biological 
years (Fig. 12A-C, Fig. 13). Biological 
year mean tail weight (Y/N,) also was 
not constant (Fig. 17B). We emphasize 
that N,/Y, and Y/N, are crude estimates 
of mean count and mean tail weight, 
respectively, and do not represent bio- 
logical year central tendency of C and w 
in the landings very well. Trends in N// 
Y. (Fig. 17A) and Y/N, (Fig. 17B) were 
cubic (sigmoid), mirroring each other 
as expected. 


Tail Weight at Which 
Half of the Biological 
Year Yield was Harvested 

The cubic trend in w50, is shown in 
Figure 18. As expected, it is similar in 
shape to that of Y/N. (Fig. 17B). How- 
ever, the two trends (Fig. 17B, Fig. 18) 
were not parallel, because the slope of 
the regression of w50; on Y/N. did not 
equal | (Fig. 19). Instead, w50, was 
1.459 times Y,/N,. Although significantly 
different from zero, the intercept of the 
regression of w50, on Y,/N, was very 
small (i.e. near the origin). 


Discussion 


It is clear that brown shrimp landings 
data should be filtered, edited, and re- 
sidual records representing outlier count 
categories removed before distributions 
of shrimp tail weight are simulated. 
The same should be (and in most cases 
have been) done before numbers of 
shrimp are estimated from landings 
within count categories, regardless of 
the algorithm used to estimate numbers 
of shrimp within count categories, 
unless the algorithms are based on 
actual sampling of size distributions 
within count categories (Ehrhardt and 
Legault, 1996). The problem of unre- 
ported landings and other limitations 
of reported landings data affect not 
only our simulations, but all other uses 
of reported landings to estimate num- 
bers of shrimp within count categories. 
These data problems cannot be rectified 
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retroactively, but should be addressed 
in the future. 

Our simulated biological year dis- 
tributions of brown shrimp tail weight 
could be biased to unknown degrees 
by many factors. This is true of all es- 
timates of numbers of shrimp derived 
from landings within count categories, 
whether at the highest possible level 
of data resolution (i.e. an individual 
shrimping trip within a statistical sub- 
area, depth zone, and month), or at lower 
levels of data resolution represented by 
various spatial-temporal aggregations 
of landings data, including ours. Our 
simulated distributions of shrimp tail 
weight should not be taken as equivalent 
to distributions of brown shrimp tail 
weight in the population of the northern 
Gulf of Mexico. However, our simulated 
distributions of w in biological year 
landings no doubt have some yet unde- 
termined relationship to actual distribu- 
tions of shrimp tail weight in the brown 
shrimp population in biological years. 
This relationship cannot be determined 
retroactively due to lack of or paucity of 
required data. Unreported landings are 
much less than reported landings, but 
our simulated distributions of shrimp 
tail weight only represent landings that 
were reported and archived. 

Despite landings data deficiencies, 
our simulated distributions of w, and 
other fishery-dependent statistics de- 
rived from them, can be useful in ex- 
amining changes in the brown shrimp 
fishery over biological years. Their 
relationships to other important fishery- 
dependent and fishery-independent 
variables could be examined in attempts 
to explain causes and effects. 

Our method could be applicable to fish- 
eries of other penaeid shrimp species for 
which landings are recorded within size 
categories expressed in C or w. It might 
also be applicable to finfish fisheries in 
which landings are reported within size 
categories expressed in number of fish 
per unit weight or in weight per fish. The 
method may also be applicable to shrimp 
landings aggregated at spatial-temporal 
levels lower (i.e. higher resolution) than 
that of an entire fishery and biological year. 

Our results suggest that brown shrimp 
were fully recruited to the fishery at 
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Figure 12.—Simulated distributions of w 
(i.e. the relationship between /’,, and w’,) 
for brown shrimp in biological years (A) 
1986, (B) 1996, and (C) 2006. 


small sizes in each biological year, then 
declined in number with w in a pattern 
similar but not identical to that of a 
negative exponential curve. In a study of 
distributions of growth rates of shrimp in 
captivity, Banks et al. (2009) examined 
effects of bin width, sample size, and 
sampling frequency on distributions 
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Figure 13.—Simulated distributions of w for brown shrimp in biological years 1986-2006, as shown with the ordinate in natural 
logarithmic scale (i.e. the relationship between In(f, ) and w’,). 
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of weight per shrimp. Interestingly, the 
shapes of their distributions of weight 
per shrimp were similar to those of our 


mean N, = 5.564x10° shrimp tails! 
No significant trend 
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Figure 14.—Simulated biological 
year total number of brown shrimp 
tails, N;, vs. coded biological year 


(T,— 1996). 
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Figure 15.—Brown shrimp yield, Y,. 
vs. coded biological year (T—1996). 
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Figure 16.—Linear relationship be- 
tween simulated total number of 
brown shrimp tails, N,, and yield, Y,. 


simulated distributions of w. Although 
we did not simulate relative distribu- 
tions of w or corresponding cumulative 
relative distributions of w, we noted that 
they could be simulated from our ap- 
proach, and they too might be of interest 
and use in shrimp stock assessments. 
Simulated biological year distribu- 
tions of w could be used to estimate 
numbers of parents and recruits, for 
purposes of determining parent-recruit 
relationships (Rothschild and Brunen- 
meister, 1984; Gracia, 1991; Ehrhardt 
and Legault, 1996; Parrack!; Nichols?). 
Numbers of parents or recruits could 
be extracted from curves representing 
distributions of w by integrating them 
over the size ranges of parents and 
recruits. However, estimates or assump- 
tions about size at maturity and growth 
patterns of males and females would 
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Figure 17.—Cubic trends in (A) 
biological year mean count, N/Y., 
vs. coded biological year (7,—-1996), 
and in (B) biological year mean 
pounds per shrimp tail, Y/N, vs. 
coded biological year (7, — 1996). 


be required, as well as estimates of 
size-specific sex ratios in the landings 
(Gracia, 1991; Ehrhardt and Legault, 
1996 Parrack', Nichols’). 

It may be possible to estimate instan- 
taneous total mortality rate (Z) from 
simulated biological year distributions 
of w by transforming them to bound- 
ed length distributions and applying 
length-based models similar to those of 
Ehrhardt and Ault (1992) (Ehrhardt®). 


°Ehrhardt, N. M. Rosentiel School of Marine 


and Atmospheric Science, Miami, FL. Personal 
commun., August 2010. 
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Figure 18.—Cubic trend in w50., 
the simulated pounds per shrimp 
tail at which half of the brown 
shrimp biological year yield, Y., was 
harvested, vs. coded biological year 
(T, — 1996). 
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Figure 19.—Linear relationship be- 
tween w50;, the simulated pounds 
per shrimp tail at which half of the 
brown shrimp biological year yield, 
Y,, was harvested, and biological 
year mean pounds per shrimp tail, 
Y/N,. 
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Alternatively, the length-based models 
used by Ehrhardt and Ault (1992) might 
be reformulated for direct application 
to biological year distributions of w for 
purposes of estimating Z (Ehrhardt®). 
Biological year distributions of w could 
also be transformed to age-frequencies 
for age-structured stock assessments. 
This would require conversion of tail 
weight to age using sex-specific growth 
curves and knowledge of size-specific 
sex ratios in the landings (Parrack, 1979; 
Rothschild and Brunenmeister, 1984; 
Gracia, 1991; Ehrhardt and Legault, 
1996; Parrack!; Nichols?). 

Simulated distributions of w of brown 
shrimp in biological year reported land- 
ings are linked, by definition and calcu- 
lations, to biological year yield. Fishing 
effort influences size-composition of the 
landings and therefore influences yield, 
although environmental variables affect- 
ing recruitment also affect yield (Cail- 
louet et al., 2008; Nance et al., 2010). 
Numbers of shrimp estimated from 
landings within count categories have 
been used in evaluating the influence 
of environmental factors on abundance, 
growth, and survival (Diop et al., 2007). 

Our method provides an alternate way 
to estimate abundance of shrimp in re- 
ported annual landings, as compared to 
algorithms used by previous investiga- 
tors. However, the relationship between 
abundance of shrimp in the landings 
and in the population remains unde- 
termined. Our simulated distributions 
of w provide examples for comparison 
with explicit or implicit assumptions 
made by previous investigators about 
the distributions of C and w. They also 
provide information of potential use in 
developing new estimators of number 
of shrimp from landings data, based on 
Statistical estimation theory and the un- 
derlying distribution of w or C. Finally, 
there may be other useful applications 
of our approach and results that we have 
not realized or anticipated. 
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Classification of Coastal Communities Reporting 
Commercial Fish Landings in the U.S. Northeast Region: 
Developing and Testing a Methodology 


SARAH L. SMITH, RICHARD B. POLLNAC, LISAL. COLBURN, and JULIA OLSON 


Introduction 


This paper introduces a method for 
classifying coastal communities for 
either sampling purposes or further 
analysis. Along the coastline from North 
Carolina to the Canadian border we find 
nearly 2,000 communities associated 
with commercial and/or recreational 
fishing. When NOAA’s National Marine 
Fisheries Service (NMFS) plans to 
implement fishery management plans, 
it is necessary to conduct (among other 
analyses) a social impact assessment 
(SIA). These SIA’s can be quite complex 
and time consuming (e.g. Pollnac et 
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al., 2006); nevertheless, they are often 
required to be submitted in a very short 
time period. In an attempt to be prepared 
to conduct SIA on short notice, all 
NMFS Regions have prepared profiles 
of a subset of the numerous coastal com- 
munities with fishing activity. These are 
called Community Profiles. This raises 
the question of how one selects the com- 
munities to be profiled. 

One hundred seventy-seven com- 
munity profiles were created and have 
been posted on the web site “Com- 
munity Profiles for the Northeast U.S. 
Fisheries” (http://www.nefsc.noaa.gov/ 
read/socialsci/community_profiles/). 
The profiles were developed as part of 
a nationwide initiative to develop com- 
munity profiles for each of the NMFS 
regions for use in Environmental Impact 
Statements (EIS). The profiles provide 
basic descriptive information, including 
a historic, demographic, cultural, and 
economic context, for understanding a 
community’s involvement in fishing and 
also furnishes a baseline from which to 
measure future change. 


ABSTRACT-— The National Marine Fish- 
eries Service is required by law to conduct 
social impact assessments of communities 
impacted by fishery management plans. 
To facilitate this process, we developed a 
technique for grouping communities based 
on common sociocultural attributes. Mul- 
tivariate data reduction techniques (e.g. 
principal component analyses, cluster an- 
alyses) were used to classify Northeast 
U.S. fishing communities based on census 
and fisheries data. The comparisons indi- 
cate that the clusters represent real group- 
ings that can be verified with the profiles. 
We then selected communities representa- 
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tive of different values on these multivari- 
ate dimensions for in-depth analysis. The 
derived clusters are then compared based 
on more detailed data from fishing commu- 
nity profiles. Ground-truthing (e.g. visiting 
the communities and collecting primary 
information) a sample of communities from 
three clusters (two overlapping geographi- 
cally) indicates that the more remote tech- 
niques are sufficient for typing the com- 
munities for further in-depth analyses. 
The in-depth analyses provide additional 
important information which we contend 
is representative of all communities within 
the cluster. 


Thus far, communities to be profiled 
have been selected on the basis of size 
and importance of fishery, types of 
fishing present, and overall knowledge 
possessed by experts working in the 
region. We posit that this technique 
is too unsystematic for this important 
endeavor, as important fishing com- 
munities could possibly be overlooked. 
SIA’s describe important implications 
of potential impacts of management ac- 
tions on fishermen and the communities 
in which they live. If SIA’s are based 
on the limited information available in 
community profiles, and if the commu- 
nities profiled are not representative of 
the communities involved in the target 
fishery, then the SIA’s produced may not 
reflect an understanding of the potential 
impact of fishery management plans 
(FMP’s). Inaccurate SIA’s can result in 
decreased fishing activity, which may 
affect household and community well- 
being and lead to social dysfunction 
within communities reliant on fishing, 
exacerbating the resistance to fisheries 
management that is evident in the North- 
east Region and elsewhere (Pollnac et 
al., 2006). 

If we could classify the large number 
of coastal communities into smaller, 
meaningful groupings, SIA data from a 
sample of communities within relevant 
subgroups would provide more accurate 
data for management decision making. 
Relevant subgroups would be those 
characterized by varying degrees of non- 
fishery and fishery attributes associated 
with participation in the target fishery or 
fisheries. Hence, the subgroups should 
be based on multivariate criteria—an 
analytic task for some form of numeri- 
cal taxonomy. 
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Figure |.—Communities profiled for the U.S. northeast community profiling project. 


Many disciplines use multivariate 
analyses for the purposes of classifica- 
tion. For example, modern biology uses 
numerical-based systematics to classify 
organisms—tools such as multiple dis- 
criminant analysis and cluster analysis. 
These techniques are not foolproof. 
First, unless all attributes of the “thing” 
to be classified are used, human decision 
making is significantly involved in the 
process. Second, a variety of techniques 
are used in numerical taxonomy (Sokal 
and Sneath, 1963), and the method se- 
lected can influence the results (Frey and 
Duek, 2007; Brusco and Kohn, 2008). 
For this reason, we felt it essential to 
test our results against several indepen- 


dent data sets, a process we refer to as 
“ground-truthing.” 


Methods 
Sample 


The attributes selected for the nu- 
merical taxonomy are derived from 
the NMFS “Social Science Data Base” 
(NMFS-SSDB) which includes com- 
mercial fisheries and U.S. Census data 
for 1,835 “ports” from North Carolina 
to the Canadian border. Those ports 
selected for community profiling are 
depicted in Figure | to demonstrate the 
geographic range of communities. By 
“ports” we mean coastal communities 


that report commercial fish landings, 
are the vessel owner port of residence, 
or homeport for permitted vessels, or 
are sites of processing, seafood/shellfish 
dealers, or recreational fishing activity. 
From the NMFS-SSDB, we selected 43 
“fishery” and 25 “social” variables for 
analysis —a total of 68 variables (Tables 
1, 2). The fishery variables selected 
were drawn from a number of variables 
characterizing fishing activity over a ten- 
year period, and included data relevant 
to quantifying fishing activity, such as 
landings by species, numbers of ves- 
sels, and numbers of vessel owners. The 
social variables used were those data 
from the 2000 United States Census that 
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could most accurately reflect changes in 
port communities that may result from or 
result in changes in fishing activity, such 
as the numbers of people employed in 
fishing related activities, the number of 
people who are self-employed, median 
household and per capita income, and 
other relevant factors. 


Data Reduction Techniques 


Principal component analysis was 
selected as the most appropriate tech- 
nique for accomplishing a reduction in 
variables because it creates a smaller 
number of new variables, grouping them 
into factors based on shared covariance. 
The 43 fishery and 25 social variables 
were reduced to fewer variables with the 
use of principal component analyses. 
The scree test (Cattell, 1966) was used 
to determine the number of components, 
resulting in four components which ac- 
count for a total of 70.4 percent of the 
variance in the data set. Components 
were rotated using the varimax tech- 
nique. The results of this analysis are 
found in Table |. Items loading highest 
on the first component (large landings, 
large vessels, sea scallops, Placopec- 
ten magellanicus; large groundfish, 
skates, Raja spp.; red crab, Geryon 
quinquedens; and monkfish, Lophius 
americanus, decreasing landings) reflect 
a fishery characterized by large vessels 
and large, but decreasing, landings of 
sea scallops, large groundfish, skates, 
red crab, and monkfish. Items load- 
ing highest on the second component 
(small vessels, many vessels, lobster, 
Homarus americanus; herring, Clupea 
harengus; and many species) indicate 
a fishery characterized by many small 
vessels, landing various species includ- 
ing lobster and herring. The third com- 
ponent reflects a fishery characterized 
by medium-sized vessels with landings 
composed principally of bluefish, Po- 
matomus saltatrix; tilefish, Lopholatilus 
chamaeleonticeps; butterfish, Peprilus 
triacanthus; mackerel, Scomber scom- 
brus; squid, Loligo pealeii, Illex illece- 
brosus; summer flounder, Paralichthys 
dentatus; scup, Stenotomus chrysops; 
and black sea bass, Centropristis striata. 
The final component refiects ports with 
changing numbers and sizes of vessels. 


Table 1.—Principal component analysis of fishery data. items in boldface type indicate highest loadings on those 
factors. 


Variable 


Value of scallops, 2003 

Landings value for home-ported vessels, 2004 
Number of large vessels (>70ft), 2004 

Average value of home-ported vessels, 1997-2003 
Value of landings at dealer reported port, 2004 
Number of large vesseis by owner city, 2003 

Total gross tonnage for home-ported vessels 
Value of large-mesh groundfish, 2003 

Value of skates, 2003 

Average landed value, 1997-2003 

Total gross tonnage for city owner vessels, 2004 
Value of red crab, 2003 

Value of monkfish, 2003 

Number of small vessels (<50ft) by owner city, 2003 
Number of small vessels by homeport, 2003 

Average number of vessels by owner city, 1997-2003 
Number of vessels by owner city, 2004 

Average number of home-ported vessels, 1997-2003 
Number of home-ported vessels, 2004 

Number of active owner city vessels, 2004 

Number of federal dealers, 2004 

Number of active home-ported vessels, 2004 
Average number of dealers, 1997-2003 

Value of lobster, 2003 

Value of herring, 2003 

Number of medium vessels (50—70ft) by owner city, 2003 
Species diversity (number of species landed), 2003 
Value of summer flounder, scup, black sea bass, 2003 
Value of butterfish, mackerel, squid, 2003 

Value of smailmesh multispecies, 2003 

Value of tilefish, 2003 

Number of medium (50-70ft) vessels by home-port 03 
Value of bluefish, 2003 

Difference in HP gross tons from 1997/98 to 2003/04 


Difference in city owner gross tons from 1997/98 to 2003/04 


Difference in HP vessels from 1997/98 to 2003/04 


Difference in number of city owner vessels from 1997/98 to 2003/04 


Value of dogfish, 2003 

Value of surf clam, ocean quahog, 2003 

Difference in dealers from 1997/98 to 2003/04 

Value of other species, 2003 

Difference in landings values for 1997/98 to 2003/04 


Difference in sum landings for HP vessels 1997/98 to 2003/04 


Percent total variance 


Component 
1 2 3 4 

0.932 0.024 0.068 0.201 
0.930 0.196 0.215 0.110 
0.932 0.184 0.219 0.084 
0.907 0.243 0.265 0.059 
0.867 0.282 0.187 0.140 
0.881 0.220 0.122 0.182 
0.852 0.345 0.326 0.154 
0.832 0.407 0.007 0.023 
0.821 0.175 0.220 0.071 
0.816 0.290 0.248 0.071 
0.789 0.376 0.185 0.302 
0.730 0.014 0.042 0.132 
0.668 0.435 0.216 —0.058 
—0.027 0.901 0.054 0.307 
0.041 0.904 0.272 0.162 
0.322 0.843 0.128 0.282 
0.350 0.825 0.105 0.346 
0.393 0.798 0.381 0.097 
0.416 0.793 0.344 0.169 
0.507 0.691 0.193 0.308 
0.487 0.657 0.071 0.020 
0.535 0.646 0.451 0.145 
0.484 0.688 0.084 -0.020 
0.087 0.575 0.012 0.091 
0.516 0.555 —0.023 —0.106 
0.502 0.525 0.281 0.282 
0.147 0.502 0.452 —0.026 
0.243 0.087 0.780 0.002 
0.193 0.103 0.710 —0.080 
0.440 0.119 0.683 0.201 
-0.093 0.070 0.648 0.431 
0.518 0.489 0.557 —0.003 
—0.007 0.107 0.488 0.147 
-0.230 0.021 -0.199 -0.776 
—0.306 —0.094 —0.339 0.656 
-0.130 —0.304 0.097 0.641 
-0.200 —-0.274 0.019 -0.622 
—0.059 0.398 0.055 0.028 
0.357 0.013 0.116 -0.006 
0.144 0.363 0.024 -0.156 
0.091 0.065 0.166 —0.025 
0.857 —0.247 0.052 —0.227 
0.928 ~0.101 —0.085 -0.242 
32.5 21.1 9.7 7.1 


Table 2 presents a principal compo- 
nent analysis of a set of variables from 
the 2000 Census. Variables selected 
can be seen in Table 2. Once again, the 
scree test was used to select number 
of components and components were 
rotated using the varimax technique. 
This resulted in three components which 
explain a total of 52.9 percent of the total 
variance in the data set. 

Component scores representing the 
position of each port on each com- 
ponent were created for each port. 
The component scores are the sum of 
the component coefficients times the 
sample standardized variables. These 
coefficients are proportional to the 
component loadings. Hence, items with 


high positive loadings contribute more 
strongly to a positive component score 
than those with low or negative load- 
ings. Nevertheless, all items contribute 
(or subtract) from the score; hence, 
items with moderately high loadings on 
more than one component (e.g. percent 
black and percent white in Table 2) will 
contribute at a moderate level, although 
differently, to the component scores as- 
sociated with each of the components. 
This type of component score provides 
the best representation of the data. 


Cluster Analysis 


Cluster analysis was then used to 
systematically group like communities 
based on these newly-created compo- 
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nent scores. As a means of combining 
the communities into relevant subgroups 
to be used for efficiently obtaining data 
for management decision making, we 
used K-means cluster analysis (Hartigan 
and Wong, 1979). The K-means proce- 
dure split the fishing communities into 
a selected number of groups by simul- 
taneously maximizing between group 
(or cluster) variation and minimizing 
within group variation. Component 
scores, which were used as input to the 
cluster analysis, are standardized, hence 
providing equal weight for each of the 
nine components used. Only cases that 
had no missing data on any of the vari- 
ables used in the principal component 
analyses are used in the cluster analysis 
(n=446). This eliminated any ports that 
did not have associated census data, 
which occurred when the port name did 
not correspond to either a geopolitically 
defined entity or a census designated 
place, bringing the number of ports 
used in the analysis from 1,835 down 
to 446. The procedure first selects the 
same number of “seeds” as the number 
of groups desired. The “seeds” selected 
are as far as possible from the center of 
all the cases. Then all cases are assigned 
to the nearest “seed,” and cases are reas- 
signed to other clusters, as needed, to 
reduce within-groups sum of squares. 

Number of clusters selected was 
based on an iterative procedure wherein 
we started at a relatively low number, 
examined the output, then increased the 
number if it was felt that, based on our 
knowledge of the ports, similar ports 
were combined. This iterative procedure 
resulted in a decision to use 40 clusters 
as the requested number. The results of 
the analysis are in Appendix I, and an 
example of selected clusters is provided 
in Table 3. 

The F-ratios across the 40 groups 
are impressive, but one must remember 
that they are an artifact of the cluster- 
ing technique which maximizes these 
values. Twelve of the clusters contain 
only one port, as illustrated by Montauk, 
N.Y., in Table 3. We believe that this is 
a valid clustering since our knowledge 
of ports included in these single port 
clusters suggests that they are unique, 
and any grouping of them with other 


Table 2.—Principal component analysis of Census data. Items in boldface type indicate highest loadings on those 


factors 


Component 


Variable 


Median household income 
High school (%) 

High school males (%) 
Poverty rate 

High school female (%) 
Unemployed (%) 

Unemployed males (%) 
Unemployed females 
Household income >200K (%) 
Share of HH income >200k 
Share of HH income retired 
Black (%) 

Males in fishing related job (%) 
Fishing related employment (%) 
Population in urban area (%) 
Females in fishing related job (%) 
Tourist housing (%) 

Hispanic (%) 

Other ethnic group (%) 

White (%) 

Two or more ethnicities 
Population 

Aggregate household income 
Asian (%) 

Male population (%) 

Percent of Total Variance 


ports would be questionable. Each of 
these single-port clusters represents a 
community with either an exceptionally 
large fishery (e.g. New Bedford, Mass.; 
Cape May, N.J.), or is a large city and 
thus the census data factors are very dif- 
ferent from the other clusters (e.g. New 
York, N.Y.; Boston, Mass.). That these 
ports appear in their own individual 
clusters indicate that they are unique 
enough to be studied on their own and 
should not be grouped with other ports. 

Note the distance for Montauk. This is 
a measure of the distance of a port from 
the center of all the cases in the cluster, 
and since there is only one, the distance 
is zero. In cluster 8, Portsmouth, N.H.., 
is closest to the center of all eight cases 
in the cluster for all seven component 
scores. Hence, this distance measure can 
be used in selecting cases from clusters 
for more intensive analysis. 

For example, one may only desire 
ports close to the center or want a 
representative sample from the cluster 
and select ports across the range of dis- 
tances. Numbers of ports in each cluster 
range from | to 123. As can be seen in 
Appendix I, many of the clusters (12) 
contain only a single case, followed 
by 7 clusters containing 2-9 cases, 2 


clusters containing 22 cases, 3 clusters 
containing 32-38 cases, and | cluster 
containing 57 cases (not all clusters are 
shown in Appendix I). 

Those clusters plotted in multidimen- 
sional space allow us to view similari- 
ties and differences on more than one 
component at a time. Figure 2 illustrates 
relative positioning of the 12 single- 
port clusters on one social component 
(population, percent in fishing related 
jobs and tourist housing) and two fishery 
components (component 2: small ves- 
sels, landing many species including 
lobster and herring and component 4: 
ports with decreasing numbers and sizes 
of vessels). A high number on fishery 
component 4 reflects rising numbers and 
sizes of vessels; hence, the name for the 
dimension— Rising. 

Figure 3 illustrates relative position- 
ing of seven multiport clusters in the 
same three-dimensional space. In this 
figure, the number following the name 
indicates cluster number as indicated in 
Appendix I. Where there are only a few 
states involved (MAME32), the states 
are abbreviated (e.g. MAME32 is cluster 
32 which includes six cases from Maine 
and Massachusetts). MIXED refers to 
too many states to abbreviate in a brief 
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0.793 0.395 0.018 

-0.766 0.172 0.413 

0.745 0.243 0.359 

0.735 0.209 0.309 

-0.732 0.088 0.444 

0.727 0.279 0.038 

id 0.659 0.277 0.029 be 

0.657 0.229 0.044 

0.624 0.302 0.104 

-0.579 0.296 0.118 

0.526 0.291 -0.247 

0.520 0.121 0.447 

0.080 0.846 0.002 

0.054 0.845 0.018 

0.156 0.599 0.271 

0.035 -0.549 0.001 

0.016 -0.475 ~0.256 

0.216 0.174 0.766 

0.276 0.135 0.745 

0.455 0.200 -0.690 

0.187 0.155 0.612 

0.078 0.095 0.570 

0.111 ~0.091 0.566 

0.230 0.280 0.451 

-0.179 -0.186 0.083 

24.072 13.358 15.469 
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title. GROUNDT refers to clusters that 
are “ground-truthed” (see below). You 
can see the ports included in cluster 8 
in Table 3. Ports included in cluster 40 
are mainly in Massachusetts with some 
from Maine, New Hampshire, and 
Rhode Island. 

Plots of clusters, such as those illus- 
trated in Figures 2 and 3, can be rotated 
to identify groups of communities that 
cluster in various selected component 
spaces, such as clusters numbered 8, 32, 
and 40. Clusters can then be examined 
by mean scores on all components, as 
in Figure 4. While communities in these 
three clusters overlap geographically 
and are quite similar on most of the 
fishery and social components, clusters 8 
and 40 are on opposite sides of the com- 
ponent mean (zero for a standardized 
variable) with regard to growth trends. 
The type of analysis presented here 
allows one to identify differences be- 
tween any subset of clusters in the data 
set, but to illustrate the process we will 
focus on these two clusters (8 and 40) 
which are used in further analyses below. 


Testing the Usefulness 
of the Cluster Analysis 


If the cluster analysis actually does 
group communities which differ on 
sociocultural and fishery variables, 
we would expect these differences 
to be manifest in other aspects of the 
community which were not measured 
as part of the original data set. To test 
this hypothesis we coded a select set 
of sociocultural variables found in the 
existing 177 community profiles, which 
were compiled from a wide range of 
available data. Eleven variables not used 
in the cluster analysis were coded, and 
percent distribution across clusters 8 and 
40 can be found in Figure 5. Despite the 
fact that there are some large differences 
between clusters 8 and 40, for example, 
presence of a fishermen’s memorial 
(50% versus 11%, respectively; Fisher’s 
Exact Test p>0.05) the small number of 
communities in each cluster (8 and 9, 
respectively) necessitates a relatively 
large difference to achieve statistical 
significance. 

It would be more revealing to exam- 
ine combinations of the sociocultural 
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Table 3.—Segment of K-Means cluster analysis output. 


Summary statistics for all cases 


Variable Between SS df Within SS df F-ratio 
FAC1FSH9 (fishery component 1) 1796.537 39 21.907 406 853.716 
FAC2FSH9 (fishery component 2) 1421.244 39 69.732 406 212.176 
FAC3FSH¢9 (fishery component 3) 1074.305 39 36.679 406 304.911 
FAC4FSHS9 (fishery component 4) 1491.341 39 100.399 406 154.635 
SOCFA‘1 (social component 1) 281.368 39 119.428 406 24.526 
SOCFA2 (social component 2) 760.698 39 88.195 406 89.790 
SOCFA3 (social component 3) 435.648 39 86.124 406 52.659 
TOTAL 7261.142 273 522.465 2842 
Cluster 7 of 40 contains 1 cases 
Members Statistics 
Case Distance Variable Minimum Mean Maximum St.Dev. 
NY, Montauk 0.00 FAC1FSH9 —4.46 —4.46 —4.46 - 
FAC2FSH9 2.60 2.60 2.60 = 
FAC3FSH9 22.73 22.73 22.73 - 
FAC4FSH9 18.56 18.56 18.56 = 
SOCFA1 0.91 0.91 0.91 = 
SOCFA2 -0.94 = 
SOCFA3 0.37 0.37 0.37 - 
Cluster 8 of 40 contains 8 cases 
Members Statistics 
Case Distance Variable Minimum Mean Maximum St.Dev. 
MA, Harwich 0.51 FAC1FSH9 —1.25 —0.64 0.62 0.61 
MA, Rockport 0.34 FAC2FSH9 3.03 4.08 5.46 0.85 
MA, Plymouth 0.52 FAC3FSH9 -0.90 0.14 0.82 0.57 
MA, Scituate 0.88 FAC4FSH9 -0.47 0.96 1.85 0.74 
ME, Kittery 0.43 SOCFA1 —0.90 0.28 0.38 0.40 
NH, Hampton 0.47 SOCFA2 -0.30 0.07 0.31 0.22 
NH, Portsmouth 0.29 SOCFA3 -0.92 -0.47 -0.20 0.24 
RI, Narragansett 0.58 
Cluster 9 of 40 contains 3 cases 
Members Statistics 
Case Distance Variable Minimum Mean Maximum St.Dev. 
ME, Stonington 0.33 FAC1FSH9 -0.85 -0.77 -0.65 0.10 
ME, Vinalhaven 0.47 FAC2FSH9 5.41 6.29 6.76 0.76 
ME, Jonesport 0.47 FAC3FSH9 —2.82 -2.47 -1.78 0.59 
FAC4FSH9 2.49 2.98 3.46 0.48 
SOCFA1 -0.33 0.33 0.73 0.57 
SOCFA2 -4.95 —4.25 -3.89 0.61 
SOCFA3 -0.10 0.13 0.36 0.23 


variables than individual items. Once 
again, we used principal component 
analysis with varimax rotation to de- 
velop scales from the profile-derived, 
sociocultural data set. Number of com- 
ponents was selected on the basis of the 
scree test. The results of the analysis are 
in Table 4. 

Table 4 indicates that the two compo- 
nents account for 43% of the variance 
in the data set. Items loading highest 
on the first component are related to 
aspects of a commercial fishing cul- 
ture, such as presence of a commercial 
fishermen’s memorial, a fishermen’s 


Table 4.—Principal component analysis of cultural and 
recreational fishing information from profiles 


Fishing Fishing 
Item Culture Recreation 
Fishermen's festival 0.667 0.258 
Blessing of fleet 0.657 -0.001 
Fishermen's memorial 0.619 —0.257 
Fishermen's assistance 0.597 0.314 
Fishermen's competition 0.553 0.107 
Fishermen's association 0.539 0.081 
Recreational fishing pier -0.090 0.718 
Fishing tournament -0.010 0.713 
Fishing education 0.361 0.487 
Percent variance 26.109 16.777 


museum, blessing of the commercial 
fleet, etc. Items loading highest on the 
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Figure 2.—Plot of single port clusters on one social and two 


fishery components. 


second component are more related to 
recreational fishing, including presence 
of a recreational fishing tournament and 
a recreational fishing pier. Presence of 
fishermen’s educational programs loads 
about the same on both components. 
Component scores, as described above, 
were calculated for each port in the 
profile data set. 

Since sample size within clusters 8 
and 40 are relatively small for statisti- 
cal analyses, we decided to cluster the 
clusters to allow comparison between 
larger groupings of ports that include 
both clusters 8 and 40. Data input were 
mean values for each of the 40 clusters 
(Appendix I) on the four fishery and 
three social component scores described 
above, and a hierarchical cluster analysis 
using median linkage and Euclidean 
distances was performed (Appendix II). 
A segment of the hierarchical tree which 
will be analyzed further is in Figure 6. 
All of the clusters found in Figure 6 can 
be found within cluster | of a K-means 
cluster analysis of the same data set 
(Appendix III). 


Fishing/Tourism 


Figure 3.—Plot of multiple port clusters on one social and 


two fishery components. 


We will now compare two clusters 
depicted in Figure 6 on the two scales 
developed from the profile data. We will 
refer to the four clusters represented by 
MASS/ME32 through MIXED38 de- 
picted at the bottom part of Figure 6 as 
Group A (n=25), and MIXED] through 
MIXED12 as Group B (n=31). Mean 
scores for Group A and Group B on the 
Fishing Culture Component are 0.297 
and —0.594, respectively (t = 4.393, df 
= 54, p<0.001), and on the Recreational 
Component they are -0.247 and 0.192, 
respectively (t= 1.581,df=54, p>0.05). 
This analysis indicates that the cluster 
analysis identified clusters that differ on 
sociocultural variables not included in 
the initial data set used for the clustering, 
providing a measure of external validity 
for the analysis. 

A final test of the usefulness of the 
clusters derived from the K-means 
cluster analysis was to “ground-truth” 
the various clusters. In contrast to the 
preceding analyses, which are based 
on secondary data (the initial database) 
and more detailed community profiles, 


which were also based on secondary 
data from publications, websites, and 
telephone inquiries as needed (see 
the community profiles), the ground- 
truthing is based on actual visits to the 
communities and interviews with com- 
munity members. 

The ground-truthing method used the 
following techniques: 


1) A photo-survey that included 
infrastructure (dock areas, fish 
processing and marketing facili- 
ties), fishing related cultural items 
(fishermen’s memorials, statues), 
and general snapshots that would 
provide an overall picture of the 
ambience of the community; 
Interviews with key informants 
concerning infrastructure and 
other points included in the pro- 
files to provide field validity 
checks; 

A brief survey that included the 
following six questions: 1) If 
you were to list five things that 
characterize [community name], 
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what would they be? 2) Would 
you say that [community name] 
is a fishing community (if not in- 
cluded in the response to the first 
question)? 3) What are three im- 
portant issues facing [community 
name] today? 4) Has [community 
name] changed over the past 5—10 
years? How? 5) Would you advise 
a young person to live in [com- 
munity name]? Why? 6) If the 
person interviewed is a fisherman, 
he or she will be asked “What’s 
it like fishing out of [community 
name]?” 


To provide a rigorous test of the 
clustering technique we selected clusters 
8 and 40 as the first two clusters to be 
compared. These two clusters overlap 
geographically and are composed of 
relatively small ports in Rhode Island, 
Massachusetts, and New Hampshire 
(Appendix I). Ground-truthed ports 
from Cluster 8 are Plymouth, Harwich, 
and Scituate, Mass., as well as Ports- 
mouth, N.H. (Sample size of surveys 
(n=89). The ports from Cluster 40 that 
were ground-truthed are Seabrook, 
N.H., and Westport, Barnstable, and 
Marshfield, Mass. (n=81). 

When ground-truthing was completed 
for the eight communities, we noted that 
communities from Cluster 8 were some- 
how “nicer.” The people in the commu- 
nities seemed to be friendlier, speaking 
of their community in a manner that 
made it seem more cohesive. These 
qualitative observations are supported 
by acontent analysis of responses made 
by community members during the 
ground-truthing exercise. While 11% 
of those interviewed in Cluster 40 said 
their communities were “spread out” 
and “composed of different parts” only 
2% of respondents from Cluster 8 made 
this observation (77 = 5.505, p<0.05). 

Additionally, a common issue in 
coastal communities is that of “gen- 
trification” —a change from being a 
fishing port to that of a desired resi- 
dential and recreational location. This 
was manifested by respondents’ com- 
plaints concerning the development of 
“condos,” “million dollar homes,” and 
an increase in “yuppies” as well as a 
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Figure 4.—Mean component values plotted for three similar clusters. 
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Figure 5.— Selected profile attributes compared across clusters 8 and 40 


“loss of character” in the port. Once 
again, Clusters 8 and 40 differed with 
respect to these responses. Forty-six 
percent of respondents from Cluster 8 
voiced these complaints in contrast to 
only 20% from Cluster 40 (y? = 13.175, 
p<0.001). These findings provide more 
external validity to the results of the 
classification methods used. 


Conclusions 


In sum, the tests of external validity 
for the cluster analyses provide support 
for the claim that the analysis actually 
did cluster communities into groupings 
that are different — different on the items 
used in the initial clustering as well as 
other variables identified by the analysis 
of the data from the community profiles 
and the ground-truthing exercise. 

We argue here that this type of 
classification of coastal communities 
is a necessary first step in providing 
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Figure 6.—Segment of hierarchial 


cluster analysis of 40 clusters from 
K-means cluster analysis 


representative information to be used 
in SIA. Community Profiles form an im- 
portant part of the information used in 
developing SIA’s, and communities to 
be profiled have thus far been selected 
on the basis of size and importance of 
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fishery, types of fishing present, and 
overall knowledge possessed by experts 
working in the region. This technique 
is too unsystematic for such an impor- 
tant endeavor. SIA’s detail important 
implications with regard to the impacts 
of management on fishermen and the 
communities in which they live. As 
noted in the introduction, the lack of 
a Statistically representative range of 
communities that may be impacted by 
proposed regulations can result in inad- 
equate SIA’s, resulting in undesirable 
effects on household and community 
well-being. All of these can exacerbate 
the types of resistance to fisheries 
management that are evident in most, 
if not all, fisheries. Using the method- 
ology described here to first select the 
communities to be profiled, as a way of 
improving the sampling process, would 
result in more representative and useful 


community profiles and, ultimately, 
improve SIA’s. 

The type of classification of coastal 
communities presented here should be 
done on a regular basis to reflect the 
rapid changes that are taking place in 
our fisheries. One of the principal com- 
ponents of the analysis of the fishery 
data reflected these changes. If regularly 
conducted, such analyses would allow 
those responsible for SIA’s to observe 
the changes in fishing communities in 
terms of their similarities and differ- 
ences, determine the factors influencing 
these changes, and use this information 
to craft more reliable and timely SIA’s 
related to specific, proposed manage- 
ment measures. 
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Appendix |.—K-Means Cluster Analysis.Distance metric is Euclidean distance, K-means splitting cases into 40 groups. Data for the following results were selected according 
to: HBOATS04> 0) AND (SOCMISDA= 0). 


Summary statistics for all cases 


Variable Between SS df Within SS 


df F-ratio 


FAC1FSH9 1796.537 39 21.907 406 853.716 
FAC2FSH9 1421.244 39 69.732 406 212.176 
FAC3FSHS9 1074.305 39 36.679 406 304.911 

FAC4FSH9 1491.341 39 100.399 406 154.635 
SOCFA1 281.368 39 119.428 406 24.526 
SOCFA2 760.698 39 88.195 406 89.790 
SOCFA3 435.648 39 86.124 406 52.659 


TOTAL 7261.142 273 522.4652 


Cluster 1 of 40 contains 57 cases 


Members 


Statistics 


Case Distance Variable 


Minimum Mean Maximum St.Dev. 

CT, Greenwich 0.58 FAC1FSH9 -0.27 -0.06 0.09 0.07 
CT, Guilford 0.21 FAC2FSH9 0.31 0.06 0.80 0.22 
CT, Madison 0.17 FAC3FSH9 0.41 -0.10 0.62 0.14 
CT, North Branford 0.22 FAC4FSH9 0.45 0.15 1.31 0.31 
MA, Aquinnah 0.52 SOCFA1 -2.98 -1.17 -0.40 0.57 
MA, West Tisbury 0.39 SOCFA2 -0.75 0.38 1.44 0.38 
MA, Georgetown 0.20 SOCFA3 -0.67 0.00 1.05 0.42 
MA, Manchester 0.41 

MA, Middleton 0.26 

MA, West Newbury 0.39 

MA, Bedford 0.17 

MA, Hopkinton 0.23 

MA, Cohasset 0.39 

MA, Dover 0.68 

MA, Norfolk 0.33 

MA, Norwood 0.26 

MA, Marion 0.46 

MA, Southborough 0.39 

MA, Sutton 0.27 

ME, Yarmouth 0.22 

NC, Ocean Island Beach 0.46 

NH, Hollis 0.25 

NH, Greenland 0.21 

NH, Hampton 0.33 

NH, New Castle 0.38 

NH, Windham 0.16 

NJ, Medford 0.12 

NJ, Avalon 0.37 

NJ, East Brunswick 0.41 

NJ, Sewaren 0.40 

NJ, Manasquan 0.27 

NJ, Monmouth 0.20 

NJ, Rumson 0.56 

NJ, Sea Bright 0.25 

NJ, Wall 0.26 

NJ, Wayne 0.24 

NY, Atlantic Beach 0.11 

NY, East Rockaway 0.23 

NY, Lido Beach 0.20 

NY, Massapequa 0.15 

NY, Seaford 0.26 

NY, Wantagh 0.20 

NY, Babylon 0.20 

NY, East Islip 0.25 

NY, Huntington Bay 0.19 

NY, Islip 0.44 

NY, Mount Sinai 0.16 

NY, Northport 0.17 

NY, Oakdale 0.23 

NY, Port Jefferson 0.20 

NY, Sayville 0.16 

NY, Southampton 0.43 

NY, Stony Brook 0.10 

NY, Armonk 0.60 

NY, Bronxville 0.89 

RI, Barrington 0.19 

Ri, East Greenwich 0.30 
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Appendix |.—(Continued). 


Cluster 3 of 40 contains 1 case 


Members 


Statistics 


Case Distance Variable Minimum Mean Maximum St.Dev. 


MA, New Bedford 0.00 FAC1FSH9 39.65 39.65 39.65 - 


FAC2FSH9 1.12 1.12 1.12 ~ 
FAC3FSH9 —0.68 -0.68 —0.68 
FAC4FSH9 6.91 6.91 6.91 _ 
SOCFA1 1.94 1.94 1.94 _ 
SOCFA2 0.10 0.10 0.10 - 


SOCFA3 


1.50 


Cluster 5 of 40 contains 1 case 


Members Statistics 


Case Distance Variable Minimum Mean Maximum 


VA, Norfolk 


0.00 FAC1FSH9 2.87 2.87 2.87 


FAC2FSH9 2.89 2.89 2.89 “ 


FAC3FSH9 3.51 3.51 3.51 - 
FAC4FSHS -15.71 -15.71 -15.71 
SOCFAI1 1.12 1.12 1.12 
SOCFA2 0.41 0.41 0.41 - 


SOCFA3 1.10 1.10 


1.10 


Cluster 6 of 40 contains 1 case 


Members 


Statistics 


Case Distance Variable Minimum Mean Maximum 


0.00 


NJ, Barnegat Light FAC1FSH9 1.27 1.27 1.27 = 


FAC2FSH9 1.67 1.67 1.67 _ 
FAC3FSH9 6.20 6.20 6.20 - 
FAC4FSH9 8.33 8.33 8.33 _ 
SOCFA1 —0.58 -0.58 -0.58 _ 
SOCFA2 -1.31 -1.31 -1.31 


SOCFA3 


Cluster 7 of 40 contains 1 case 


Members Statistics 


Case Distance Variable Minimum 


Mean Maximum 


NY, Montauk 


0.00 FAC1FSH9 


FAC2FSH9 2.60 2.60 2.60 - 
FAC3FSH9 22.73 22.73 22.73 _ 
FAC4FSH9 18.56 18.56 18.56 — 
SOCFA1 0.91 0.91 0.91 - 
SOCFA2 -0.94 - 


SOCFA3 


0.37 


Cluster 8 of 40 contains 8 cases 


Members Statistics 


Case Distance Variable Minimum 


Mean Maximum 


MA, Harwich FAC1FSH9 0.64 0.62 0.61 


MA, Rockport 0.34 FAC2FSH9 3.03 4.08 5.46 0.85 


= MA, Plymouth 0.52 FAC3FSH9 -0.90 -0.14 0.82 0.57 
3 MA, Scituate 0.88 FAC4FSH9 ~0.47 0.96 1.85 0.74 
% ME, Kittery 0.43 SOCFA1 -0.90 -0.28 0.38 0.40 
NH, Hampton 0.47 SOCFA2 -0.30 0.07 0.31 0.22 
NH, Portsmouth 0.29 SOCFA3 -0.92 0.47 -0.20 0.24 
Ri, Narragansett 0.58 
continued 
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Appendix |.—(Continued). 


Cluster 10 of 40 contains 1 case 


Members 


Case 


Distance 


Statistics 


Variable 


Minimum 


Mean 


Maximum 


NY, New York 


0.00 FAC1FSH9 0.61 0.61 0.61 
FAC2FSH9 4.31 4.31 4.31 - 
FAC3FSH9 ~0.83 - 
FAC4FSH9 -4.77 -4.77 -4.77 
SOCFA2 -4.73 -4.73 -4.73 - 
SOCFA3 15.44 15.44 15.44 


Cluster 11 of 40 contains 32 cases 


Members 


Case 


Distance 


Statistics 


Variable 


Minimum 


Mean 


Maximum 


DE, Frederica 
DE, Milford 
DE, Millsboro 
MA, Onset 
MD, Cambridge 
MD, Crisfield 

MD, Willards 

MD, Berlin 

MD, Snow Hill 

ME, Eastport 

NC, Aurora 

NC, Belhaven 

NC, Gloucester 
NC, Marshallberg 
NC, Morehead City 
NC, Newport 

NC, Swan Quarter 
NC, Wilmington Beach 
NC, Bayboro 

NC, Vandemere 
NJ, Millville 

NJ, Keansburg 

NJ, Neptune City 
NY, Mastic Beach 
Ri, East Providence 
Rl, Woonsocket 
VA, Melfa 

VA, Onancock 

VA, Hallwood 

VA, Exmore 

VA, Nassawadox 
VA, Portsmouth 


FAC1FSHS9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 
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Appendix |.—(Continued). 


Cluster 12 of 40 contains 38 cases 


Statistics 


Distance Variable Minimum Mean Maximum 


, Lewes 0.34 FAC1FSH9 0.55 
Brewster 0.34 FAC2FSH9 0.66 1.58 
Dennis 0.39 FAC3FSH9 0.39 
Eastham 0.31 FAC4FSHS 0.97 2.31 
South Dennis 0.23 SOCFA1 0.21 0.64 
. Yarmouth 0.42 SOCFA2 K -0.09 0.82 
Vineyard Haven 0.32 SOCFA3 0.04 
. Essex 0.25 
Newburyport 0.71 
Salisbury 0.40 
Swampscott 0.38 
Nantucket 0.38 
Kingston 0.24 
Middleboro 0.26 
Ocean Bluff 0.33 
. Falmouth 0.37 
Scarborough 0.47 
South Portland 0.37 
Hancock 0.43 
Buxton 0.34 
Kittery 0.34 
Ogunauit 0.46 
Saco 0.32 
, Wells 0.34 
Newington 
Dover 
Middietown 
, Beach Haven 
Forked River 
, Manahawkin 
Point Pleasant 
J, Toms River 
Tuckerton 
Waretown 
NY, Oceanside 
Ri, Charlestown 
VA, Wachapreague 
VA, Poquoson 


Cluster 14 of 40 contains 9 cases 


Statistics 


Case Distance Variable Minimum Mean Maximum 


MA, Gosnold 0.70 FAC1FSH9 -0.07 -0.01 
MD, Smith Island 0.49 FAC2FSHS —0.26 0.20 
ME, Cranberry Isles 0.43 FAC3FSHS -0.27 -0.09 
ME, Matinicus 0.43 FAC4FSH9 0.14 0.15 0.68 
ME, North Haven 0.51 SOCFA1 -1.05 0.09 1.89 
ME, Roque Bluffs 0.64 SOCFA2 -6.51 -5.05 -3.37 
NC, Smyrna 0.43 SOCFA3 0.04 0.62 1.84 
VA, Saxis 0.82 

VA, Tangier 0.51 


Cluster 15 of 40 contains 1 case 


Statistics 


Case Distance Variable Minimum Mean Maximum 


MA, Gloucester 0.00 FAC1FSH9 1.53 1.53 1.53 
FAC2FSH9 25.40 25.40 25.40 
FAC3FSH9 —2.45 —2.45 -2.45 
FAC4FSH9 1.57 1.57 1.57 
SOCFA1 -0.03 -0.03 
SOCFA2 -0.01 
SOCFA3 -0.30 


continued 


Marine Fisheries Review 


a 
Members 
0.24 
0.50 
0.33 
0.49 
0.44 
0.41 
: 0.33 
vu 
St.Dev. 
0.05 
0.14 
0.08 
0.25 
: 0.96 
1.06 
0.63 
Members 
St.Dev. 
i 
52 


Appendix |.—(Continued). 


Cluster 16 of 40 contains 1 case 


Statistics 


Case Distance 


Variable Minimum Mean Maximum St.Dev. 


NJ, Cape May 0.00 FAC1FSH9 8.09 8.09 


8.09 - 
FAC2FSH9 1.28 1.28 1.28 - 
FAC3FSH9 6.75 6.75 6.75 — 
FAC4FSH9 3.01 3.01 3.01 - 
SOCFA1 0.40 0.40 0.40 - 
SOCFA2 0.07 0.07 0.07 ~ 
SOCFA3 


Cluster 17 of 40 contains 1 case 


Members Statistics 


Case Distance Variable Minimum Mean Maximum St.Dev. 


MA, Chatham 0.00 FAC1FSH9 -2.58 -2.58 -2.58 
FAC2FSH9 12.65 12.65 12.65 - 
FAC3FSH9 1.28 1.28 1.28 - 
FAC4FSH9 0.82 0.82 0.82 - 
SOCFA1 -0.17 -0.17 
SOCFA2 -0.73 -0.73 -0.73 ~ 


SOCFA3 


-0.37 -0.37 


-0.37 


Cluster 18 of 40 contains 1 case 


Members Statistics 


Case Distance Variable Minimum Mean Maximum St.Dev. 


ME, Portland 0.00 FAC1FSH9 5.16 5.16 5.16 - 
FAC2FSH9 11.07 11.07 11.07 - 
FAC3FSH9 -0.70 -0.70 -0.70 - 
FAC4FSH9 -14.62 -14.62 -14.62 - 
SOCFA1 0.09 0.09 0.09 - 
SOCFA2 0.27 0.27 0.27 - 


SOCFA3 


0.12 


0.12 


Cluster 19 of 40 contains 22 cases 


Members Statistics 


Case Distance 


Variable Minimum Mean 


Maximum 


CT, Bridgeport 0.66 FAC1FSH9 0.15 0.01 0.36 


0.11 
CT, Norwalk 0.40 FAC2FSH9 -0.27 0.02 0.48 0.20 
CT, Stamford 0.45 FAC3FSH9 0.23 0.02 0.30 0.15 
CT, New Haven 0.62 FAC4FSH9 -1.18 0.31 0.33 0.44 
DE, Wilmington 0.60 SOCFA1 0.48 0.67 2.13 0.80 
MA, Lynn 0.34 SOCFA2 0.04 0.52 1.19 0.26 
MA, Framingham 0.49 SOCFA3 0.73 1.72 3.67 0.81 
MA, Randolph 0.47 
MA, Revere 0.49 
MA, Worcester 0.17 
NJ, Ventnor City 0.30 
NJ, Jersey City 0.77 
NJ, Long Branch 0.24 
NJ, Clifton 0.23 
NY, Baldwin 0.44 
NY, Glen Cove 0.27 
NY, Inwood 0.36 
NY, Staten Isiand 0.39 
NY, Bay Shore 0.26 
PA, Philadelphia 0.60 
RI, Providence 0.73 


VA, Richmond 0.47 


Cluster 20 of 40 contains 1 case 


Members Statistics 


Case 


Distance Variable Minimum Mean Maximum 


ME, Harpswell 0.00 FAC1FSH9 -1.75 -1.75 -1.75 


FAC2FSHS9 5.45 5.45 5.45 
FAC3FSH9 -2.57 -257 -2.57 - 
FAC4FSH9 9.44 9.44 9.44 - 
SOCFA1 -0.57 -0.57 -0.57 
SOCFA2 -1.58 -1.58 -1.58 
SOCFA3 -0.21 0.21 - 


continued 


4 

Members 

| 

: 

St.Dev. 

{ 
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Cluster 23 of 40 contains 22 cases 


Case Distance 


Statistics 


Variable Minimum 


Maximum 


DE, Bowers 
MA, Chilmark 
ME, Brookline 
, Brooksville 
, Castine 
Franklin 
. Sorrento 
Sullivan 
Tremont 
Isle au Haut 
, St. George 
, Bremen 
Bristol 
Southport 
=, Georgetown 
Columbia 
Jonesboro 
Harkers Island 
, Ocracoke 
, Sneads Ferry 
Orient 
Onley 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


Cluster 25 of 40 contains 35 cases 


Case Distance 


Statistics 


Variable Minimum Mean 


Maximum 


CT, Branford 0.26 
CT, East Lyme 
CT, Groton 

CT, Mystic 

CT, Noank 

MA, Danvers 
MA, Ipswich 

MA, Methuen 

MA, Nahant 

MA, Salem 

MA, Saugus 

MA, Quincy 

MA, Weymouth 
MA, Duxbury 
MA, Hingham 
MA, Hull 

MA, Pembroke 
ME, Cape Elizabeth 
ME, Bath 

ME, Eliot 

ME, Kennebunk 
ME, York 

ME, York Harbor 
NJ, Atlantic City 
NJ, Beimar 

NJ, Brielle 

NY, Island Park 
NY, Point Lookout 
NY, East Hampton 
NY. East Quogue 
NY, West Islip 

Ri, Warwick 

Ri, Jamestown 
RI, Cranston 

Ri, Westerly 


FAC1FSH9 -0.47 0.13 
FAC2FSH9 

FAC3FSH9 

FAC4FSHS 

SOCFA1 

SOCFA2 

SOCFA3 
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Members 
Mean — St.Dev. 

0.15 0.33 -0.10 0.04 0.11 

0.58 0.28 0.13 1.38 0.48 

0.17 ~0.43 -0.09 0.09 0.11 i 

0.18 0.06 1.07 0.37 
0.35 -0.79 0.87 0.46 

0.39 -2.75 -1.70 0.56 

0.37 -1.07 0.53 0.11 0.30 

: 0.28 

0.44 

j 0.24 

: 0.34 

0.49 

0.42 

0.33 

0.45 

0.36 

0.17 

: 0.52 

0.22 

0.49 

0.42 

0.36 

Members 

0.14 0.14 
1.45 0.35 

0.66 0.29 

0.18 0.43 

0.33 0.49 
0.82 0.27 

0.70 0.38 
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Cluster 28 of 40 contains 1 case 


Members 


Case 


Distance 


Statistics 


Variable Minimum 


Maximum 


MA, Boston 


0.00 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


Cluster 31 of 40 contains 1 case 


Members 


Distance 


Variable Minimum 


Maximum 


0.00 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


Cluster 32 of 40 contains 6 cases 


Members 


Case 


Distance 


Statistics 


Variable Minimum 


Maximum 


MA, Orleans 

MA, Truro 

MA, Wellfleet 

ME, Bar Harbor 

ME, Southwest Harbor 
ME, Boothbay Harbor 


0.52 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


Cluster 35 of 40 contains 8 cases 


Members 


Case 


Distance 


Variable Minimum 


Maximum 


DE, Leipsic 

MA, Buzzards Bay 
ME, Gorham 

ME, Machias 

NC, Elizabeth City 
NH, Durham 

NJ, Wildwood 

RI, Kingston 


FAC1FSH9 
FAC2FSH9 

FAC3FSH9 

FAC4FSH9 

SOCFA1 

SOCFA2 

SOCFA3 


0.00 


Cluster 36 of 40 contains 6 cases 


Members 


Case 


Distance 


Statistics 


Variable Minimum 


Mean Maximum 


CT, Stonington 
MA, Falmouth 
NJ, Sea Isle City 
NY, Mattituck 

RI, Little Compton 
RI, Tiverton 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


continued 


5.47 5.47 5.47 ; 
0.25 0.25 0.25 
-7.88 -7.88 -7.88 
0.32 0.32 0.32 
0.35 0.35 0.35 
2.75 2.75 2.75 
Statistics 
VA, Newport 5.30 5.30 5.30 
-2.58 -2.58 -2.58 
1.52 1.52 1.52 = 
5.13 5.13 5.13 - 5 
0.61 0.61 0.61 = 
0.34 0.34 0.34 
0.83 0.83 0.83 
; -1.09 ~0.55 0.18 0.32 
ae 0.65 1.59 2.35 3.01 0.55 e 
0.53 ~0.70 0.29 0.41 
0.32 0.71 0.24 1.27 0.71 
0.42 0.24 0.66 1.49 0.62 3 
0.44 -1.02 0.48 0.65 
-1.64 1.11 ~0.68 0.38 
Statistics 
0.26 ~0.04 = 0.04 
0.33 -0.09 0.16 0.15 a 
0.29 -0.07 0.17 0.11 
0.43 0.16 0.30 
0.56 2.00 2.96 0.52 eo 
0.12 1.21 2.84 0.76 
3 0.21 0.50 
0.64 
0.68 ~0.38 0.87 0.51 
0.40 0.26 0.65 1.16 0.35 o 
0.50 0.83 1.48 2.46 0.72 e 
0.35 -1.44 ~0.25 0.74 0.79 “ 
4 0.63 ~0.69 0.11 0.72 0.57 oe 
0.40 -0.71 0.06 0.81 0.52 
-1.06 0.77 0.22 
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Cluster 38 of 40 contains 4 cases 


Case 


Distance 


Statistics 


Variable Minimum Mean 


Maximum 


MA, Sandwich 
NC, Beaufort 
VA, Chincote 
VA, Virginia 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 


Cluster 40 of 40 contains 9 cases 


Members 


Case 


Distance 


Statistics 


Variable Minimum Mean 


Maximum 


MA, Barnstable 
MA, Westport 
MA, Beverly 
MA, Marblehead 
MA, Newburyport 

, Marshfield 
ME, Kennebunkport 
NH, Rye 
NH, Seabrook 


FAC1FSH9 
FAC2FSH9 
FAC3FSH9 
FAC4FSH9 
SOCFA1 
SOCFA2 
SOCFA3 
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Members 
0.48 0.63 0.33 0.42 
0.81 0.56 1.84 2.94 1.00 
0.50 0.85 1.62 2.35 0.61 4 
0.72 0.68 1.20 2.41 0.82 hoes 
-0.51 0.33 1.04 0.81 
-1.20 0.28 0.14 0.62 
~0.66 -0.05 1.11 0.79 
0.29 ~0.84 0.56 ~0.40 0.14 
0.45 2.05 2.78 3.44 0.48 *¢ 
0.32 -0.01 0.34 0.96 0.37 
0.52 -1.52 -0.71 0.17 0.62 
0.35 -1.40 -0.34 0.24 0.49 
0.37 -0.26 0.32 0.76 0.37 
0.35 
0.36 
56 


Appendix Il. 


of K-means 40 Clusters. Distance metric is Euclidean distance. Median 
linkage method. Single port clusters use port name rather than cluster number. 


Cluster 
containing 


Cluster 


containing 


Were joined 
at distance 


No. of members 
in new cluster 


MIXED25 
MIXED1 
MIXED1 
MIXED1 
MIXED36 
MIXED36 
MIXED36 
GROUNDT40 
GROUNDT8 
MIXED35 
GROUNDT8 
GROUNDT8 
GROUNDT8 
GROUNDT8 
GROUNDT8 
CTNJ39 
MAINE33 
GROUNDT8 
MIXED14 
MIXED14 
MIXED14 
MIXED14 
MIXED14 
NY22 
HARPSWELL 
NY22 
NEWPORT VA 
BARNG LIGHT 
NCNY13 
NCNY13 
NCNY13 
NCNY13 
PORTLAND 
NCNY13 
NCNY13 
NEW YORK 
NEW YORK 
NEW YORK 


NEW BEDFORD 


MIXED21 
MIXED25 
MIXED12 
MIXED23 
MIXED1 
MIXED11 
MIXED19 
MAME32 
GROUNDT40 
MIXED36 
MIXED38 
MIXED35 
MIXED26 
MIXED24 
MENC37 
NJRI4 
MAINES 
MAME34 
MAINE2 
GROUNDT8 
NCNJ29 
MANJ30 
CTNJ39 
MIXED14 
MAINE33 
BOSTON 
CAPE MAY 
NEWPORT VA 
BARNG LIGHT 
NY22 
HARPSWELL 
MAINE27 
NORFOLK 
CHATHAM 
PORTLAND 
NCNY13 
GLOUCESTER 
MONTAUK 
NEW YORK 


0.411 
0.378 
0.437 
0.599 
0.624 
0.631 
0.639 
0.674 
0.665 
0.773 
0.797 


0.779 13 
0.850 14 
0.773 15 
0.945 16 
1.023 2 
1.196 2 
1.256 17 
1.353 2 
1.294 19 
1.482 20 
3 1.516 21 
1.772 23 
1.593 24 
2.329 3 
2.351 25 
2.835 2 
2.578 3 
2.541 4 
2.918 29 
2.660 32 
2.725 33 
3.657 2 
3.703 34 
4.775 36 
4.930 37 
7.041 38 
11.146 39 
| 
| 
{ 
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NEW BEDFORD 
GLOUCESTER 
CAPE MAY 
NEWPORT VA 
BARNG LIGHT 


MIXED14 
MAINE2 
GROUNDTS 
GROUNDT40 
MAME32 
MIXED38 
MIXED19 
MIXED11 
MIXED23 
MIXED12 
MIXED21 
MIXED25 
MIXED1 
MIXED36 
MIXED35 
MIXED26 
MIXED24 
MENC37 
MAME34 
NY22 
BOSTON 
MAINES 
MAINE33 
HARPSWELL 
MAINE27 
CHATHAM 
PORTLAND 
NORFOLK 
NEW YORK 
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NJRI4 ~------------------------------------------- 
CTNJ39 
MANJ30 
= 
+-- 
‘ +----------------------- 
+----- 
MONTAUK 
; 
| > 


Appendix Iil.—K-means clustering of K-means 40 clusters. K-means splitting 40 cases into 10 groups (single port clusters use port name rather than cluster number). 


Summary statistics for all cases 


TOTAL 


4492.010 


511.949 


Variable Between SS df Within SS df F-ratio 

LARGE 1590.686 a 86.966 30 60.970 

SMALL 784.370 9 88.921 30 29.403 

MEDIUM 674.282 9 61.448 30 36.577 
: RISING 1039.884 9 178.613 30 19.407 
se POVERTY 28.208 9 19.093 30 4.925 
fe URBAN 141.997 9 40.185 30 11.779 
E ETHNIC 232.583 9 36.724 


Cluster 1 of 10 contains 21 cases 


Members 


Case 


Statistics 


Distance 


Variable 


Minimum 


Mean 


Maximum 

MIXED1 0.80 LARGE -0.64 

MAINE2 1.34 SMALL -1.33 1.09 4.13 1.42 
Rr GROUNDT8 1.30 MEDIUM -0.99 0.20 2.11 0.78 

MIXED11 0.62 RISING -2.92 -0.19 3.78 1.56 

MIXED12 0.62 POVERTY -1.17 0.44 2.06 0.84 

MIXED19 0.81 URBAN -2.74 -0.19 1.21 0.96 

MIXED21 0.52 ETHNIC -1.13 0.06 4.37 1.22 

NY22 2.03 

MIXED23 0.75 

MIXED24 1.14 

MIXED25 0.48 

MIXED26 0.91 

NCNJ29 1.58 

MANJ30 1.77 

MAME32 0.74 

MAME34 1.66 

MIXED35 1.02 

MIXED36 0.62 

MENC37 0.99 

MIXED38 0.82 


GROUNDT40 


Cluster 2 of 10 contains 3 cases 


Members 


Case 


Statistics 


Distance 


Variable 


Minimum 


Maximum 


NORFOLK 
PORTLAND 2.18 
BOSTON 2.06 


LARGE 


SMALL 2.89 
MEDIUM -0.70 
RISING -15.71 
POVERTY 0.09 
URBAN 0.27 


ETHNIC 


Cluster 3 of 10 contains 1 case 


Members 


Case 


Distance 


Variable 


Minimum 


Maximum 


NEW BEDFORD 0.00 


LARGE 39.65 
SMALL Li 
MEDIUM —0.68 
RISING 6.91 
POVERTY 1.94 
URBAN 0.10 
ETHNIC 1.50 


continued 


a 
0.83 
1.84 3.29 5.16 1.70 
6.48 11.07 4.18 ae 
1.02 3.51 2.21 oe 
-12.74 -7.88 4.24 
0.51 1.12 0.54 
0.34 0.41 0.07 i 
0.12 1.24 2.75 1.44 
Statistics 
1.12 1.12 - 
-0.68 
6.91 6.91 
‘ 1.94 1.94 
0.10 0.10 - 
1.50 1.50 - : 
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Cluster 4 of 10 contains 6 cases 


Members Statistics 


Case Distance Variable Minimum Mean Maximum 


NJRI4 LARGE 
BARNG LIGHT SMALL 
NCNY13 4 MEDIUM 
CAPE MAY RISING 
NEWPORT VA POVERTY 
CTNJ39 URBAN 
ETHNIC 


Cluster 5 of 10 contains 1 case 


Members 


Case Distance Variable Minimum Maximum 


CHATHAM 0.00 LARGE -2.58 
SMALL 
MEDIUM 
RISING 
POVERTY 
URBAN 
ETHNIC 


Cluster 6 of 10 contains 1 case 


Members 


Case Distance Variable Minimum Maximum 


MONTAUK 0.00 LARGE —4.46 
SMALL 2.60 
MEDIUM 22.73 
RISING 18.56 
POVERTY 0.91 
URBAN 
ETHNIC 0.37 


Cluster 7 of 10 contains 1 case 


Members Statistics 


Case Distance Variable Minimum Mean Maximum 


NEW YORK 0.00 LARGE 
SMALL 
MEDIUM 
RISING 
POVERTY 
URBAN 
ETHNIC 


Cluster 8 of 10 contains 3 cases 


Members Statistics 


Case Distance Variable Minimum Maximum 


MAINE 1.11 LARGE 

HARPSWELL 1.71 SMALL 

MAINE33 0.97 MEDIUM 
RISING 
POVERTY 
URBAN 
ETHNIC 


continued 
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-1.16 2.31 8.09 3.61 
-2.58 0.42 1.67 1.54 
1.52 5.88 10.15 2.80 
3 -2.18 2.69 8.33 3.76 ce 
: -0.58 0.07 0.61 0.41 
0.21 0.34 0.74 
0.08 0.91 0.67 
: Statistics 
St.Dev. 
-2.58 -2.58 A 
a 12.65 12.65 
0.82 0.82 
-0.17 -0.17 ~ 
-0.37 -0.37 = 
Statistics 
2.60 2.60 
22.73 22.73 
18.56 18.56 
0.91 0.91 - 
~0.94 = 
0.37 0.37 
St.Dev. 
0.61 0.61 0.61 - 
4.31 4.31 4.31 
-4.77 -4.77 -4.77 
5 15.44 15.44 15.44 = 
StDev. 
-1.75 -1.36 0.52 
: 4.50 5.41 6.29 0.90 
-2.57 -2.26 -1.73 0.46 
2.98 5.64 9.44 3.38 
~0.57 0.33 0.46 
-5.89 -3.91 -1.58 2.18 
0.19 0.64 0.43 
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Cluster 9 of 10 contains 2 cases 


Case 


Variable 


Minimum 


MIXED14 
MAINE27 


LARGE 
SMALL 
MEDIUM 
RISING 
POVERTY 
URBAN 
ETHNIC 


-9.43 
0.62 


Cluster 10 of 10 contains 1 case 


Case 


Variable 


Minimum 


GLOUCESTER 


LARGE 
SMALL 
MEDIUM 
RISING 
POVERTY 
URBAN 
ETHNIC 


1.53 
25.40 


73(2) 


Members Statistics 
0.97 -0.07 -0.05 ~0.03 0.03 
0.97 -0.09 0.02 0.14 0.16 ; 
3 -0.09 -0.04 0.02 0.08 os 
-0.37 -0.11 0.15 0.37 
‘ -2.35 -1.13 0.09 1.73 2 
-7.24 -5.05 3.10 
1.04 1.45 0.59 
2 
4 Members Statistics 4 
Distance Mean Maximum St.Dev. 3 
25.40 25.40 
-2.45 -2.45 -2.45 
1.57 1.57 1.57 
-0.03 -0.03 -0.03 
-0.01 -0.01 -0.01 
-0.30 -0.30 -0.30 
t 
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The Marine Fisheries Review publishes re- 
view articles, original research reports, signifi- 
cant progress reports, technical notes, and 
news articles on fisheries science, engineering, 
and economics, commercial and recreational 
fisheries, marine mammal studies, aquaculture, 
and U.S. and foreign fisheries developments. 
Emphasis, however, is on in-depth review ar- 
ticles and practical or applied aspects of ma- 
rine fisheries rather than pure research. 

Preferred paper length ranges from 4 to 12 
printed pages (about 10-40 manuscript pages), 
although shorter and longer papers are some- 
times accepted. Papers are normally printed 
within 4-6 months of acceptance. Publication 
is hastened when manuscripts conform to the 
following recommended guidelines. 


The Manuscript 


Submission of a manuscript to the Marine 
Fisheries Review implies that the manuscript 
is the author’s own work, has not been sub- 
mitted for publication elsewhere, and is ready 
for publication as submitted. Commerce De- 
partment personnel should submit papers 
under a completed NOAA Form 25-700. 

Manuscripts must be typed double-spaced 
throughout and submitted with two duplicate 
copies. The complete manuscript normally 
includes a title page, a short abstract, text, 
literature citations, tables, figure legends, foot- 
notes, and the figures. The title page should 
carry the title and the name, department, insti- 
tution or other affiliation, and complete ad- 
dress (plus current address if different) of the 
author(s). Manuscript pages should be num- 
bered and have |-inch margins on all sides. 
Running heads are not used. An “Acknowl- 
edgments” section, if needed, may be placed 
at the end of the text. Use of appendices is 
discouraged. 


Abstract and Headings 


Keep titles, headings, subheadings, and the 
abstract short and clear. Because abstracts are 
circulated by abstracting agencies, it is impor- 
tant that they represent the research clearly 
and concisely. Headings within each section 
must be short, reflect a logical sequence, and 
follow the rules of multiple subdivision (i.e. 
there can be no subdivision without at least 
two items). 


Style 


The entire text should be intelligible to inter- 
disciplinary readers; therefore, all acronyms, 


Editorial Guidelines for 


the Marine Fisheries Review 


abbreviations, and technical terms should be 
spelled out the first time they are mentioned. 
The scientific names of species must be writ- 
ten out the first time they are mentioned; sub- 
sequent mention of scientific names may be 
abbreviated. Follow the U.S. Government 
Printing Office Style Manual (1984 ed.) and 
the CBE Style Manual (Sth ed.) for editorial 
style, and the most current issue of the Ameri- 
can Fisheries Society’s Common and Scien- 
tific Names of Fishes from the United States 
and Canada for fish nomenclature. Only jour- 
nal titles, scientific names (genera and spe- 
cies), and vessel names should be italicized. 
Dates should be written as follows: 11 Nov. 
199i. Measurements should be expressed in 
metric units, e.g. metric tons as t; other equiva- 
lent units may also be listed in parenthesis. 
Common abbreviations and symbols such as 
mm, m, g, ml, mg, and °C (without periods) 
may be used with numerals. The numeral one 
(1) should be typed as a one, not as a lower- 
case el (1). Write out the numbers zero through 
nine unless they form part of measurement 
units (e.g. nine fish but 9 mm). 


Footnotes 


Footnotes should not be embedded within 
the text document. They must be numbered 
with Arabic numerals and typed on a separate 
sheet of paper. Footnote all personal commu- 
nications, listing the name, affiliation, and ad- 
dress of the communicator and date of com- 
munication. Unpublished data and unpub- 
lished manuscripts should include the title, 
author, pagination of the manuscript or re- 
port, and the address where it is on file. Au- 
thors are advised to avoid references to non- 
standard (gray) literature, such as internal, 
project, processed, or administrative reports, 
wherever possible. Where these references are 
used, please include whether they are avail- 
able from NTIS (National Technical Informa- 
tion Service) or from some other public de- 
pository. 

Literature Cited 


Title the list of references “Literature Cited” 
and include only published works or those 
actually in press. Citations must contain the 
complete title of the work, inclusive pagina- 
tion, full journal title, and the year, month, 
volume, and issue numbers of the publication. 
Unpublished reports or manuscripts and per- 
sonal communications must be footnoted. 

Citations should be double-spaced and listed 


alphabetically by the senior author’s surname 
and initials. Coauthors should be listed by ini- 
tials and surname. Where two or more cita- 
tions have the same author(s), list them chro- 
nologically; where both author and year match 
on two or more, use lowercase alphabet to 
distinguish them (1979a, 1979b, 1979c, etc.). 
Authors must double-check all literature cited; 
they alone are responsible for its accuracy. 


Tables 


Tables should be printed separately and 
double-spaced. Tables should not be exces- 
sive in size and must be cited in numerical 
order in the text. Headings should be short but 
sufficent to allow the table to be intelligible on 
its own. All unusual symbols must be ex- 
plained in the table heading. Other incidental 
comments may be footnoted with Arabic nu- 
merals. Because tables are typeset, they need 
only be submitted typed and formatted, with 
double-spaced legends. Zeros should precede 
all decimal points for values less than one. 
Table headings and format should be consis- 
tent; do not use vertical rules. 


Figures 


Figures include line illustrations and photo- 
graphs (or slides) and must be cited in numeri- 
cal order in the text. Figures are to be labeled 
with author’s name and number of figure. Use 
Times Roman font (upper and lowercase let- 
ters) to label within figures. Avoid vertical let- 
tering except for y-axis labels. Zeros should 
precede all decimal points for values less than 
one. Figures should be submitted as both la- 
ser-printed copies and computer software files. 
Figure legends should explain all symbols and 
abbreviations and should be double-spaced on 
a separate page at the end of the manuscript. 
Consider column and page sizes when design- 
ing figures. Please note that we do not print 
graphics in color. 


Finally 


First-rate, professional papers are neat, ac- 
curate, and complete. Authors should proof- 
read the manuscript for typographical errors 
and double-check its contents and appearance 
before submission. Mail the manuscript to: 


Editor, Marine Fisheries Review 

Scientific Publications Office 

National Marine Fisheries Service, NOAA 
7600 Sand Point Way N.E. 

Seattle, WA 98115 
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